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ABSTRACT 


/ 

Five  base-line  calibrations  of  the  Motorola  Mini-Ranger 
III  (MRS  III)  shortwave  positioning  system  were  performed 
across  Puget  Sound,  Washington,  in  January  1984.  Two  MRS 
III  codes  were  calibrated  over  measured  distances  of  1061.2, 
2417.5,  4083.1,  7016.8,  and  9861.1  meters.  Reference  lengths 
were  obtained  to  submeter  accuracy  using  a  Tellurometer , 

Model  CA-1000 . 

Results  of  the  data  analysis  suggest  that  an  optimum 
base-line  length  of  about  2,500  meters  provides  a  character¬ 
istic  calibration  curve  that  adequately  depicts  all  equip¬ 
ment  configurations  for  a  typical  MRS  III.  The  calibration 
curve  was  best  characterized  by  high-gain  antenna  data  for 
two  Mini-Ranger  codes.  The  high-gain  antenna  also  exhibited 
the  least  sample  deviation  of  the  three  antenna  types  cali¬ 
brated.  Range  average  was  found  to  be  insignificant.  Re¬ 
gression  curves  applied  to  the  data  showed  no  change  in 
quality  of  fit  above  a  power  of  three  and  variances  of  0.5 
meters  indicated  a  fair  to  good  fit  of  the  data. 
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I.  INTRODUCTION 

A.  BACKGROUND  INFORMATION 

Electronic  positioning  of  hydrographic  survey  vessels  is 
a  common  practice  in  most  automated  hydrographic  agencies. 

The  National  Ocean  Service  (NOS)  is  such  an  organization  and 
is  responsible  for  charting  the  United  States  territorial 
waters.  The  liability  that  accompanies  this  responsibility 
has  prompted  NOS  to  analyze  the  quality  of  its  positioning 
data.  Specifically,  NOS  would  like  to  apply  f unctionally- 
derived  range  error  correctors  to  range  rates  obtained  from 
superhigh  frequency  (SHF)  positioning  systems  in  real  time. 
SHF  positioning  systems,  such  as  Motorola's  Mini-Ranger  III 
(MRS  III),  are  used  for  inshore  and  harbor  surveys  of  scales 
greater  than  1:20,000.  Portability,  accuracy,  and  automation 
are  some  of  the  reasons  hydrographers  prefer  them. 

The  systems  generally  read  out  directly,  in  meters,  the 
distance  or  range  of  the  vessel.  However,  these  range  rates 
are  not  error  free  due  to  weather  changes,  attenuation  dis¬ 
tances,  changing  components  in  a  given  configuration,  and 
other  variables.  The  errors  are  determined  through  base¬ 
line  calibrations  where  observed  ranges  are  compared  against 
a  known  distance  and  error  curves  are  generated  that  display 
range  error  as  a  function  of  signal  strength.  Range  correc¬ 
tors  are  then  extracted  from  these  curves  and  applied  to  the 
received  range  rates  in  near-real  time. 
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In  January  1984,  calibrations  were  performed  over  five 


measured  base  lines  that  extended  across  Puget  Sound, 
Washington.  Three  different  antenna- types  and  three  differ 
ent  Range  Average  settings  were  sampled  for  each  base  line. 
The  results  were  used  to  explain  the  observed  relationship 
between  range  error  and  signal  strength  over  a  variety  of 
conditions  and  system  configurations. 

1.  SHF  Positioning  Systems 


SHF  positioning  systems  have  played  a  major  role  in 
the  development  of  automated  hydrography  in  NOS.  In  general, 
electronic  methods  have  replaced  antiquated  manual  methods 
for  positioning  on  scales  as  large  as  1:10,000.  These  micro- 
wave  "black  boxes"  were  for  years  rejected  by  stubborn  "moss- 
back"-  hydrographers  who  refused  to  accept  the  accuracy 


claims . 


Today's  electronic  systems  have  freed  the  hydro- 


grapher  of  the  tedium  of  manually  acquiring  and  processing 
data  so  the  hydrographer  may  concentrate  on  planning,  error 
analysis,  data  interpretation,  and  other  more  important 
functions.  MRS  III  has  been  a  mainstay  in  west  coast  sur¬ 
vey  operations  in  NOS  since  the  early  1970's. 

Though  most  of  the  problems  encountered  with  MRS  III 
have  been  solved  in  the  course  of  the  agency's  experience 
with  the  system,  the  problem  of  systematic  range  errors  and 
their  necessary  correction  has  not  been  fully  resolved. 


Additionally,  a  new  automated  hydrographic  data  acquisitic 


W 


system  called  SDS  III  is  slated  to  replace  the  present  HYDRO¬ 
PLOT/HYDROLOG  system  in  1987.  The  SDS  III  development  group 
is  concerned  about  present  SHF  calibration  methods  and  would 
like  to  devise  a  procedure  for  accurately  determining  range 
correctors  over  a  variety  of  conditions  and  applying  them  in 
near-real  time. 

2 .  Theory  of  Operation 

The  theory  of  SHF  positioning  systems  is  an  exchange 
of  pulsed  signals,  where  radar-like  signals  are  transmitted 
from  the  vessel  to  an  active  shore  station.  Upon  interroga¬ 
tion,  the  shore  station  returns  a  signal  to  the  vessel. 
Distance  is  directly  related  to  the  round-trip  travel  time 
of  the  signal  which  is  based  on  a  fixed  speed  of  propagation 
of  the  electromagnetic  waves.  Atmospheric  effects  will  vary 
this  speed  of  propagation  but  corrections  can  be  applied 
given  the  proper  meteorological  data.  [Ref.  l:pp.  45-52] 

In  the  specialized  case  of  the  MRS  III,  the  system 
uses  one  or  more  reference  stations  placed  over  points  of 
known  position  on  shore.  Each  reference  station  can  trans¬ 
mit  different  pulse  forms  (called  codes).  These  codes  can 
be  set  in  the  reference  stations  and  then  selected  for  in¬ 
terrogation  by  setting  a  dial  on  the  range  console.  The 
reference  stations  are  actively  interrogated  by  a  receiver- 
transmitter  assembly  ( RT  unit)  located  aboard  the  vessel. 

The  range  console  (also  aboard  the  vessel)  contains  a  stable 


counter  clock  that  measures  the  travel  time  and  converts  it 


to  a  distance.  The  range  console  also  performs  some  internal 
smoothing  and  screening  functions  on  the  received  signals. 

The  resultant  range  rates  are  displayed  at  the  front  of  the 
console.  The  system  comes  with  three  different  types  of  an¬ 
tennae  whose  purposes  are  discussed  in  greater  detail  in 
Chapter  II.  The  MRS  III  operates  on  a  frequency  of  5.4  to 
5.6  GHz  (versus  9  GHz  for  most  SHF  systems).  This  lower  fre¬ 
quency  helps  reduce  the  possibility  of  interference  from 
other  radars  in  the  working  area.  The  system  has  been  used 
by  NOS  for  hydrographic  surveys  ranging  from  scales  of 
1:20,000  to  as  large  as  1:5,000,  when  used  in  a  hybrid  mode 
such  as  range-azimuth. 

3 .  Field  Uses  of  SHF  Systems 

NOS  has  been  using  MRS  Ill's  for  over  a  decade  in  a 
variety  of  survey  applications.  Most  positioning  schemes 
use  the  MRS  III  in  either  a  range-range  or  range-azimuth 
geometry.  Field  edit  positions,  velocity  casts,  and  bottom 
samples  may  also  be  positioned  this  way. 

SHF  systems  are  very  portable  and  are  well-suited 
for  inshore  hydrography  and  shoreline  development  requiring 
a  shallow-draft  vessel.  Also,  the  positions  of  control 


stations  in  a  project  area  can  be  verified  by  performing  a 
system  check  calibration  of  the  MRS  III  in  the  field.  Sys¬ 
tem  checks  are  discussed  in  greater  detail  later  in  this 
chapter . 


The  systems  are  also  used  in  NOS  to  check  lane  counts 
of  medium-frequency  positioning  systems  such  as  RAYDIST.  In 
NOS,  MRS  III  range  rates  are  recorded  at  periodic  intervals 
to  explain  or  pinpoint  where  RAYDIST  lane  jumps  (lane  losses 
or  gains)  have  occurred.  Once  lane  jumps  have  been  deter¬ 
mined,  whole  lane  correctors  can  be  added  to  existing  partial 
correctors  of  the  RAYDIST  to  allow  the  hydrographer  to  con¬ 
tinue  survey  operations.  This  procedure  has  saved  consider¬ 
able  amounts  of  hydrographic  data  by  relieving  the 
hydrographer  of  the  need  to  return  to  a  calibration  or  ini¬ 
tialization  area. 

SHF  systems  are  practical  for  positioning  a  vessel 
using  range-azimuth  geometry.  This  geometry  uses  a  shore- 
based  observer  with  a  theodolite  positioned  next  to  the  SHF 
shore  station.  The  observer  radios  observed  azimuths  to  the 
survey  launch.  This  setup  is  less  probable  with  a  medium- 
frequency  system  like  RAYDIST  because  the  observer's  pres¬ 
ence  in  the  ground  plane  of  the  antenna  affects  the  outgoing 
signal  significantly.  The  popularity  of  range-azimuth  geom¬ 
etry  has  prompted  the  recent  development  of  automatic  and 
semi-automatic  total  stations  that  consist  of  an  SHF  unit 
and  electronic  theodolite.  Accuracy  is  increased  by  using 
an  azimuth  vice  a  second  range  rate  in  the  position 
computation . 


There  are  several  options  that  can  greatly  enhance 
the  operation  of  the  MRS  III  and  help  reduce  the  effect  of 


systematic  range  errors: 

a.  A  signal  strength  indicator  displays  the  strength  of 
each  incoming  signal.  Signal  strengths  can  be  recorded 
in  real  time  and  used  to  identify  position  blunders  or 
busts.  They  can  also  be  used  to  determine  range  error 
correctors  in  near  real  time. 

b.  A  range  averaging  option  enables  the  user  to  select 
the  number  of  range  rates  to  be  averaged  before  dis¬ 
playing.  This  option  was  originally  installed  to 
provide  more  accuracy  by  increased  averaging  of  the 
data. 

c.  A  multi-user  capability  offers  16  different  codes. 

The  extra  codes  give  the  hydrographer  flexibility  to 
work  in  areas  where  other  Mini-Rangers  are  operating. 
Though  not  directly  related  to  range  errors,  anomalous 
range  rates  can  sometimes  be  attributed  to  other  MRS 
III  systems  located  in  the  survey. 

d.  A  voltage  cut-off  switch  powers  the  unit  down  when  a 
minimum  of  22  volts  is  not  exceeded.  It  was  found 
that  voltages  below  22  volts  caused  significant  range 
errors. 

e.  The  Space  Diversity  option  solves  the  major  problem 
of  range  holes  or  null  zones  by  using  two  RT  units. 

The  nulling  phenomenon  is  prevented  by  mounting  two 
RT  unit  antennae  at  different  heights  on  the  survey 
vessel  with  a  switch  to  select  between  the  antennae, 
until  an  acceptable  signal  is  received. 


B.  SYSTEM  ERRORS 

1 .  Systematic  Range  Errors 

MRS  III,  though  an  accurate  system  (Appendix  A)  to 
use,  still  contains  "some  idiosyncracies  that  have  remained 
largely  unknown  and  which  can  lead  to  serious  range  error 
problems."  [Ref.  2:p.  5]  This  statement  was  made  in  1972. 


Today  many  of  the  error  sources  have  been  resolved  but  others 
are  as  mysterious  as  ever.  Operation  of  the  system  can  be 
affected  by  range,  meteorological  conditions,  range  holes  or 
null  zones,  antenna  types,  range  averaging,  power  supply 
voltage,  warm-up  time,  and  even  the  RT  unit  cable  length. 

Only  errors  due  to  range,  range  averaging,  and  antenna  types 
were  considered  in  evaluating  the  January  1984  calibration 
data . 

Changes  in  weather  within  a  survey  project  area  can 
cause  drift  in  the  MRS  III.  Temperature  changes  directly 
affect  a  unit's  circuitry--specif ically  its  magnetron. 

"Very  few  weather  phenomena  directly  or  significantly  affect 
the  range  measurement  of  the  Mini-Ranger  system  because  of 
the  high  frequency  of  operation--5  GHz.  However,  meteorolo¬ 
gical  conditions  do  contribute  to  variations  of  the  effec¬ 
tive  path  length  (apparent  range)  and,  more  severely, 
atmospheric  absorption  (attenuation)."  [Ref.  3: Section  3.1.2] 

The  actual  propagation  speed  of  the  radio  frequency 
(RF)  waves  is  a  function  of  wavelength,  air  temperature, 
barometric  pressure,  and  vapor  pressure.  The  counter  clock 
inside  the  range  console  uses  a  fixed,  average  value  for  the 
propagation  speed.  Deviations  from  the  fixed  value  are 
determined  during  a  base-line  calibration  performed  in  the 
field  by  the  user. 

Phenomena  of  "high  particle  density  such  as  rain, 
mist,  fog,  or  photochemical  smog"  are  also  responsible  for 


attenuation.  A  formula  (Eqn.  1.1)  exists  for  calculating  the 
attenuation  due  to  rain  [Ref.  4:p.  11]  where  A  is  attenuation 
in  db,  r  is  rainfall  rate  in  mm/hour,  and  a  &  b  are  tabulated 
frequency  functions. 

A  =  ar*3  (1.1) 


Signal  attenuation  reduces  the  range  of  operation  and 
lowers  signal  strengths,  so  calibrations  should  be  performed 
as  near  to  the  working  area  as  possible  to  minimize  errors 
due  to  weather.  Thus,  any  range  errors  found  during  the 
base-line  calibration  will  be  a  function  of  weather  condi¬ 
tions  in  the  project  area.  Table  I  shows  attenuation  esti¬ 
mates  derived  by  using  Equation  1.1. 


TABLE  I 


ATTENUATION  ESTIMATES  FOR  A  3 -CM  WAVE  LENGTH 


Precipitation 
Rate  (mm/hr ) 

6.0 

15.2 
18.7 
22.6 

34.3 
43.1 


Attenuation 
( db/km ) 

0.192 

0.613 

0.728 

0.801 

1.280 

1.640 


Systematic  range  errors  are  the  focus  for  this  thesis 
Their  presence  is  attributed  to  the  fundamental  radar  equa¬ 
tion  (Eqn.  1.2),  where  P'  is  the  intensity  of  radiated  power, 
P  is  the  peak  power,  G  is  the  antenna  gain,  and  R  is  the 


distance  from  the  source. 


The  effects  of  the  radar  equation  are  shown  by  plotting  range 
error  versus  distance  in  kilometers.  The  same  curve  form  is 


also  found  in  plots  of  range  error  versus  signal  strength. 

In  the  latter  case,  the  errors  result  from  artificially- 
attenuated  signals  that  simulate  greater  distances.  Thus, 
NOS  base-line  calibrations  of  the  MRS  III  are  related  to  the 
radar  equation.  In  a  typical  NOS  calibration  curve  (Fig. 
1.1),  the  range  error  curve  becomes  exponential  below  a  cer¬ 
tain  critical  signal  strength  threshold.  At  high  signal 
strengths  the  curve  exhibits  a  near  linear  form.  Casey  re¬ 
fers  to  the  threshold  value  separating  the  two  curve  types 
as  the  Critical  Strength  Threshold  (CST).  NOS  calls  the  CST 
the  Minimum  Acceptable  Signal  Strength  (MASS)  and  NOS  re¬ 
jects  data  below  this  value.  Casey  has  found  that  above  the 
CST  "the  errors  are  distributed  with  a  Gaussian  or  normal 
distribution;  below  the  CST  the  distribution  is  log  normal." 
[Ref.  2:p.  9] 

The  important  point  is  if  the  curve  were  considered 
"well-defined"  below  the  MASS,  data  acquired  with  signal 
strengths  below  this  threshold  might  not  be  rejected.  This 
amounts  to  a  large  savings  in  time  and  money  when  viewed 
from  a  (NOS)  fleet  aspect.  A  well-defined  curve  is  one 
whose  data  exhibit  low  standard  deviations. 
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Range  errors  are  due  principally  to  pulse  shape  deg 
radation  of  the  received  signals.  To  determine  the  amount 
of  degradation  and  the  effect  on  the  received  range  rate,  a 
signal  strength  integer  value  is  used  that  defines,  through 
calibration,  what  range  corrector  (the  algebraic  inverse  of 
range  error)  should  be  applied  to  the  range  rate. 

2 .  Signal  Strength  as  an  Indicator 

NOS  MRS  Ill's  are  equipped  with  the  optional  signal 
strength  indicator  that  displays  a  "relative  number  output 
that  is  directly  proportional  to  the  amplitude  of  the  re¬ 
ceived  signal  from  the  reference  station."  [Ref.  3:Section 
5.2.3]  Signal  strengths  are  measured  because  tests  have 
shown  that  range  error  is  dependent  on  signal  strength  in  a 
repeatable  manner.  Signal  strengths  are  significant  because 
the  counter  clock,  which  measures  the  elapsed  time,  must 
recognize  the  leading  edge  of  the  received  pulse  because 
".  .  .  the  leading  edge  shape  of  the  incoming  pulse  varies 
with  the  amplitude  so  that  a  low  amplitude  (read  signal 
strength)  leads  to  a  systematic  timing  (read  range)  error." 
[Ref.  4:p.  9]  This  relationship  between  elapsed  counter  time 
and  amplitude  is  best  shown  by  diagram  (Fig.  1.2). 

The  factors  affecting  MRS  III  signal  strength  that 


are  discussed  in  this  thesis  are  range  averaging,  base-line 
lengths,  and  three  types  of  reference  station  antennae. 


C.  NOS  CALIBRATION  PROCEDURES 


1.  Base-Line  Calibrations 


Base-line  calibrations  are  required  by  NOS  (1)  at 
the  beginning  and  end  of  a  survey  project,  (2)  at  intervals 
not  exceeding  2  months  during  long  projects,  (3)  whenever 
an  RF  component  is  changed  in  either  an  RT  unit  or  tran¬ 
sponder,  and  (4)  whenever  the  difference  between  prior  base¬ 
line  correctors  and  daily  system  check  correctors  exceeds 
specified  rejection  limits-  [Ref.  5:p.  M-2] 

A  description  of  the  base-line  calibration  procedure 
is  taken  from  the  PMC  OPORDER,  Appendix  M: 

The  base-line  calibration  is  performed  by  introducing 
signal  attenuation  with  a  variable  attenuator  over  a 
known  distance.  The  observed  range  readings  are  compared 
to  the  known  distance  at  various  signal  strengths  to 
determine  range  errors.  A  graph  is  constructed 
which  typically  shows  a  nearly-linear  portion  followed  by 
a  sharp  rise  at  low  signal  strength  values.  A  mean  error 
is  determined  for  signal  strength  values  corresponding 
to  the  nearly-linear  portion  of  the  graph  and  its  al¬ 
gebraic  inverse  is  the  corrector  to  be  applied  to  the 
observed  range  readings.  The  signal  strength  value 
corresponding  to  the  sharp  upswing  portion  of  the  curve 
is  the  minimum  value  allowed  for  hydrographic  survey 
operations.  [Ref.  5:pp.  M-2  &  M-3 ] 

For  a  base-line  calibration,  ten  range  rates  are  sampled  for 
each  signal  strength  to  obtain  the  mean  range.  Base-line 
lengths  are  measured  prior  to  calibration  to  submeter  accu¬ 
racy  using  an  NOS-approved  laser  or  microwave  electronic 
distance  measuring  instrument. 

NOS  presently  requires  the  length  of  the  base  line 


to  be  between  700  and  3,000  meters.  The  optimum  length  has 
yet  to  be  determined  and  is  discussed  in  the  thesis. 


Though  the  correctors  determined  from  a  base-line 
calibration  are  most  often  those  applied  to  the  final  posi 


tion  data,  NOS  does  require  that  a  closing  or  bracketing 
calibration  be  performed  at  the  completion  of  a  survey. 
Additionally,  daily  critical  or  noncritical  system  checks 
are  made  to  monitor  any  drift  in  the  Mini-Rangers.  These 
system  checks  test  for  instrument  drift  and  verify  or  dis¬ 
pute  the  most  recent  base-line  correctors. 

2 .  Critical  System  Checks 

A  critical  system  check  is  defined  as  one  in  which 
an  independent  method  of  position  determination  is  employed 
to  check  range  errors.  Noncritical  system  checks  are  much 
less  time  consuming  and  are  performed  in  conjunction  with  a 
recent  critical  system  check.  This  approach  was  devised  by 
NOS  to  minimize  the  calibration  time  and  maximize  the  time 
spent  acquiring  hydrographic  data.  Though  noncritical  sys¬ 
tem  checks  require  less  time,  they  do  not  provide  an  inde¬ 
pendent  check  fix  to  verify  the  first  fix.  Whenever  the 
difference  between  a  daily  system  check  and  base-line  cali¬ 
bration  corrector  exceeds  the  rejection  limits  (Table  II) 
and  cannot  be  resolved  after  a  careful  check  is  made,  a  new 
base-line  calibration  is  required  to  continue  hydrographic 
survey  operations.  [Ref.  5:pp.  M-5  &  M-6] 

The  following  are  methods  currently  approved  for 


critical  system  checks. 


TABLE  II 


NOS  CALIBRATION  REJECTION  LIMITS 


Survey  Scale 


Rejection  Limit  (m) 


1:  10,000  5 

1:  20,000  and  smaller  10 


a.  Fixed-Point  Method 


The  fixed-point  method  is  the  critical  system 
check  method  preferred  by  NOS.  It  involves  bringing  the 
launch  alongside  a  geodetically-positioned  pile,  wharf,  rock, 
or  other  point  feature  and  comparing  the  observed  range  rates 
with  precomputed  range  rates.  It  is  operationally  and  com¬ 
putationally  simplest  of  all  the  methods  and  is  not  restrict¬ 
ed  by  visibility  in  the  work  area.  However,  the  method  is 
not  always  feasible  as  frequently  a  point  feature  is  not 
available  in  the  work  area.  Additionally,  this  method  is 
restricted  to  launches  since  deep-drafted  survey  ships  can¬ 
not  afford  the  grounding  risk  of  the  shoal  waters  surround¬ 
ing  such  features. 

b.  Navigational  Range  and  Cut-Off  Angle  Method 
The  navigational  range  and  a  cut-off  angle 
method  can  be  employed  when  as  few  as  two  geodetic  control 
stations  are  in  the  area.  To  construct  a  navigational 
range,  one  geodetic  control  station  is  used  as  a  front  (or 
rear)  range  marker.  A  second  station,  or  erected  signal, 
is  used  as  the  other  range  marker. 


the  above  method  (b).  Though  this  may  require  considerable 
set-up  time,  the  method  is  quick  and  convenient  once  the 
ranges  are  in  place.  To  calibrate,  the  launch  steers  one 
range  and  crosses  the  other.  The  observed  range  rates  are 
compared  with  precomputed  rates  by  subtraction. 

e.  Three-Point  Sextant  Resection  Method 

The  three-point  sextant  resection  method  has 
been  used  extensively  in  NOS  pending  the  development  of  more 
efficient  methods.  It  requires  a  minimum  of  four  visual 
targets  and  three  sextant  observers.  Simultaneous  to  the 
reading  of  the  MRS  III  range  rates,  horizontal  sextant  an¬ 
gles  are  observed  to  the  four  signals.  The  observed  sextant 
angles  are  input  into  a  calibration  program.  Though  compu¬ 
tations  are  performed  by  an  onboard  computer,  they  are  still 
tedious  and  more  involved  than  the  simpler  subtractions  used 
in  some  other  methods. 

f.  Theodolite  Intersection 

The  theodolite  intersection  method  requires 
three  observers,  each  with  a  theodolite  to  turn  angles  from 
an  initial  station  to  the  survey  launch.  The  recorded  azi¬ 
muths  define  a  fix  and  check  fix  that  are  input  into  a  cali¬ 
bration  program.  Although  it  is  one  of  the  more  accurate 
methods,  it  is  rarely  used  by  NOS  personnel  due  to  the 
equipment  and  labor  demands.  The  ensuring  computations  are 
only  slightly  faster  than  for  the  resection  method. 


Noncritical  system  checks  are  generally  quicker  than 
critical  system  checks  but  do  not  include  an  independent 
position  determination.  Additionally,  they  cannot  check  geo¬ 
detic  control  in  the  working  area.  NOS  has  approved  their 
use  because  the  MRS  III  has  been  used  in  west  coast  survey 
operations  for  over  a  decade  and  has  developed  into  a  relia¬ 
ble  system  that  usually  does  not  develop  flaws  which  cannot 
be  detected  using  noncritical  system  checks  in  conjunction 
with  periodic  critical  system  checks.  This  fact,  coupled 
with  the  policy  of  using  base-line  correctors  in  lieu  of 
separate  correctors  determined  daily,  has  been  the  driving 
force  behind  reducing  the  requirements  for  twice-daily  criti¬ 
cal  system  checks.  [Ref.  5:p.  S-2) 
a.  Three-Range  Method 

The  three-range  method  checks  the  operation  of 
three  separate  shore  stations  by  comparing  the  three  possible 
fixes  generated  by  the  three  range  rates.  A  major  constraint 
is  that  the  vessel  be  positioned  to  obtain  a  near  60-degree 
angle  of  intersection  for  each  of  the  three  lines  of  position 
(LOPs)  to  obtain  consistently  strong  fixes  for  each  pair. 
Unlike  the  other  two  noncritical  system  check  methods  that 
utilize  ready-made  calibration  forms,  the  three-range  method 
requires  separate  fix  computations  and  an  onboard  computer. 
The  key  to  using  this  method  is  realizing  that  it  is  sensi¬ 
tive  to  the  intersection  angles  of  the  LOPs. 


b.  Base-Line  Crossing  Method 

The  base-line  crossing  method  compares  a  sum  of 
two  observed  range  rates  to  the  known  base-line  length  con¬ 
necting  two  transponder  stations.  Thus,  the  launch  crosses 
between  two  shore  stations  and  a  mark  is  called  when  a  shore 
observer  radios  that  the  crossing  launch  is  on  range  with  the 
two  shore  stations.  To  stabilize  the  update  rate  of  the  MRS 
III,  the  launch  steers  a  range  arc  of  cne  of  the  stations. 

The  ensuing  computations  are  subtractive  and  simple  to 
perform. 

c.  Launch  to  Launch  Method 

The  launch  to  launch  method  is  very  quick  because 
many  launches  can  be  checked  at  once.  Using  a  critically- 
checked  launch,  designated  the  guide  launch,  other  launches 
come  alongside  with  their  bows  pointed  at  the  transponder  to 
be  checked.  They  do  not  have  to  be  tied  but  they  should  be 
aligned  along  their  beams.  MRS  III  rates  are  simultaneously 
recorded  and  the  guide  launch  correctors  applied.  Non-guide 
launches  then  apply  their  respective  base-line  correctors 
and  compare  to  the  rejection  limits  table  (Table  II).  All 
system  check  methods,  both  critical  and  noncritical,  require 
the  hydrographer  to  first  apply  the  most  recent  base-line 
correctors  to  the  observed  rates  before  comparing  them 
against  their  known  values. 


D.  FOCUS  OF  RESEARCH 


1 .  The  Future  of  SHF  Units 

SHF  positioning  systems  have  been  available  for  over 
a  decade.  They  may  eventually  be  replaced  by  futuristic 
satellite  systems  such  as  the  NAVSTAR  Global  Positioning 
System  (GPS).  The  need  to  locate  transponders  on  shore  for 
launch  operations  could  be  eliminated  along  with  the  some¬ 
times  complex  geodetic  networks  that  support  a  survey.  With 
only  6  of  the  scheduled  18  satellites  in  orbit,  testing 
has  already  begun  to  determine  GPS  accuracies.  The  entire 
system  is  expected  to  be  operational  in  1987.  Other  advan¬ 
tages  to  this  powerful  system  are  its  continuous  passive 
output,  multiple  user  capabilities,  and  world-wide  coverage. 
I-t  will  also  free  the  hydrographer  from  any  constraints  due 
to  fix  geometry  or  line-of-sight .  Studies  have  been  con¬ 
ducted  by  the  Naval  Oceanographic  Office  (NAVOCEANO)  in  Bay 
St.  Louis,  Mississippi,  to  determine  GPS's  suitability  for 
dynamic  positioning  of  a  survey  launch.  Thus  far,  GPS  accu¬ 
racies  of  10  meters  have  been  observed  in  a  dynamic  mode. 

[ Ref .  6 : p.  4  ] 

A  cost  study  recently  done  by  NAVOCEANO  [Ref.  7:p. 
84]  states  a  savings  of  nearly  four  million  dollars  a  year 
for  two  Navy  hydrographic  survey  ships  would  be  a  result  of 
using  GPS.  It  appears  likely  that  the  GPS  system  will  make 
dinosaurs  of  not  only  SHF  positioning  systems,  but  medium- 
frequency  systems  as  well.  If  GPS  receivers  were  made 


affordable  even  to  fishermen,  it  would  replace  LORAN-C  and 
Omega.  However,  there  are  several  reasons  why  the  present 
systems  will  still  be  used: 


GPS  accuracies  for  hydrography  (a  dynamic  mode)  are  not 
yet  verified.  This  fact  alone  suggests  that  SHF  systems 
will  continue  to  be  used  in  NOS  for  large-scale 
hydrography . 

GPS  may  not  be  able  to  keep  pace  with  the  short  fix  in¬ 
terval  used  to  obtain  hydrographic  positions  on  large- 
scale  surveys.  Also,  additional  computer  storage  may 
be  necessary  to  accomodate  the  large  amounts  of  GPS 
data . 

Access  to  the  precise  (P)  code  may  be  restricted  to 
military  operations.  Since  NOS  does  not  directly  sup¬ 
port  such  operations,  it  may  not  be  granted  access  to 
this  code.  Access  to  the  P  code  also  generates  con¬ 
siderable  security  problems  regarding  its  use  by  NOS 
personnel.  If  access  to  the  P  code  is  limited  there 
are  other  techniques,  such  as  interferometry  and 
differential  GPS,  that  will  improve  the  accuracy  of 
data  obtained  with  the  coarse  acquisition  (CA)  code. 

The  initial  cost  for  a  GPS  receiver  is  prohibitive  to 
the  average  mariner  and  may  also  be  too  expensive  for 
NOS.  "Initial  estimates  of  costs  are  $50,000  for 
each  system  (excluding  the  cost  of  the  receivers) 
and  $15,000  to  $30,000  one-time  cost  for  tailored 
software."  [Ref.  7:p.  83]  Using  this  estimated,  one 
can  imagine  the  initial  cost  to  NOS  for  nearly  50 
operating  systems. 

2 .  Improving  MRS  III  Calibration  Methods 

The  above  statements  suggest  that  MRS  III  and  similar 
systems  will  still  be  used  during  the  evolution  and  testing 
of  SDS  III.  Base-line  calibrations,  to  determine  systematic 
range  errors,  will  still  be  necessary.  The  present  base-line 
calibration  methods  can  be  automated  and  enhanced  to  acquire 
more  data  in  less  time  without  manual  manipulation  by  the  ob¬ 
server.  One  such  practice  is  to  record  calibration  data 
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directly  into  a  portable  computer.  The  same  program  used  to 
acquire  the  data,  for  example  LOTUS  1-2-3,  can  be  used  to 
automatically  average  the  correctors  and  graph  them.  NOAA 
Ship  DAVIDSON  is  presently  using  this  procedure  for  base¬ 
line  calibrations  of  its  Mini-Ranger  Falcon  484,  the  MRS  III 
successor.  Possible  refinements  to  the  method  include  a 
mathematical,  nonsubjective ,  determination  of  the  MASS  and 
applying  curve— fitting  routines  to  obtain  a  functional  form 
of  the  calibration  curve. 

A  further  enhancement  would  be  to  tabulate  the  cor¬ 
rectors  into  a  corrector  file  such  as  the  Sound  Velocity, 
Tide,  Zone,  Parameter,  and  Station  Tables  that  exist  pres¬ 
ently  in  the  Bathymetric  Swath  Survey  System  (BSSS). 

Another  possibility  is  to  incorporate  the  functional  form 
(equation)  of  the  calibration  curve  directly  into  the 
start-up  routine  of  the  SURVEY  or  HYDROPLOT  acquisition  pro¬ 
grams  by  entering  the  equation's  coefficients.  This  could 
easily  be  implemented  by  prompting  the  operator  for  the 
equation's  coefficients.  The  creation  of  a  more  automated 
method  for  base-line  calibrations  could  provide  a  corrector 
curve  and/or  tape  for  direct  application  to  the  acquired 
navigation  data.  Automation  of  the  current  base-line  cali¬ 
bration  procedures  and  the  ensuing  computations  that  yield 
the  calibration  curve  would  enable  the  operator  to  sample 
more  calibration  data  in  the  same  amount  of  time. 


3.  Present  Calibration  Problems 


Some  issues  need  to  be  addressed  regarding  present 
NOS  methods  of  determining  range  error  correctors  for  cor¬ 
responding  signal  strengths: 

The  assumption  of  a  best-fit  linear  curve  being  most 
representative  of  an  MRS  Ill's  range  error-signal 
strength  relationship  is  not  necessarily  the  most 
accurate  one  considering  the  variety  of  conditions 
and  equipment  configurations  over  which  a  system  can 
operate.  Another  choice  is  to  use  a  separate  cor¬ 
rector  for  each  signal  strength  in  the  calibration 
curve.  This  approach  assumes  that  a  functional  re¬ 
lationship  exists  between  range  error  and  signal 
strength  and  can  be  defined  through  curve-fitting 
procedures.  This  thesis  discusses  what  degree 
polynomial  best  fits  the  January  1984  data  after 
curve-fitting  techniques  were  tried. 

The  optimum  base-line  length  for  a  base-line  cali¬ 
bration  has  been  revised  several  times  by  NOS. 

As  new  field  phenomena  are  realized  the  recommended 
calibration  length  changes  [Ref.  8].  When  base 
lines  as  short  as  100  meters  are  used,  corrector 
anomalies  occur  at  the  higher  signal  strengths. 

The  converse  situation  of  observing  range  anomalies 
at  low  strengths  for  longer  base  lines  limits  the 
base-line  length  on  the  high  end.  The  cause  of 
the  anomalies  it  close-in  ranges  is  "the  measured 
pulse  amplitude  riding  on  top  of  the  previously 
decaying  pulse."  [Ref.  8]  Anomalies  occur  at 
ranges  less  than  700  meters  where  long  distances 
are  being  simulated  with  a  variable  attenuator. 

The  dilemma  is  if  too  short  a  base-line  length 
is  selected,  the  resultant  signal  strength  versus 
range  corrector  curve  may  not  become  exponential  at 
low  signal  strengths  and  thus  MASS  cannot  be  determined. 
On  the  other  hand,  if  too  long  a  base-line  length 
is  selected,  no  correctors  will  be  obtained  at 
higher  signal  strengths  and  correctors  for  field 
data  must  be  extrapolated  from  lower  signal  strength 
correctors  (Ref.  9:p.  2]. 

At  present  in  NOS,  MASS  is  a  subjective  estimate  of 
the  point  of  inflection  (or  signal  strength)  where 
the  curve  becomes  exponential.  A  mathematical 
definition  of  MASS  that  can  be  used  in  future  cali¬ 
bration  software  is  desirable.  Although  the  aim  is 


to  fit  a  continuous  function  to  the  calibration  curve, 
the  curve  is  really  not  a  discrete  one  in  that  non¬ 
integer  signal  strengths  do  not  exist.  In  the  case  of 
the  newer  Falcon  system,  every  signal  strength  may  not 
be  calibrated  due  to  the  expanded  signal  strength 
scale.  The  curve  is  thus  empirically  derived  from 
those  strengths  that  were  calibrated.  The  MRS  III 
signal  strength  scale  is  from  0  to  about  25,  while 
the  Falcon  scale  range  is  from  0  to  100. 

Present  base-line  calibration  procedures  are  not  re¬ 
sponsive  to  observed  anomalies  in  the  calibration 
curve  because  a  single  mean  corrector  is  selected 
for  the  entire  curve.  By  using  a  functional  form  of 
the  entire  curve,  anomalies  could  either  be  rejected 
through  an  outlier  statistical  test  or  become  incor¬ 
porated  into  the  correction  process.  The  presence 
of  these  anomalies  is  confirmed  by  data  obtained  for 
this  thesis.  Anomalies  are  shown  in  the  range-error 
versus  signal-strength  plots  in  Appendix  B. 

This  research  should  aid  NOS  in  developing  a  calibration 


procedure  that  realistically  predicts  the  range  corrector 
for  an  observed  range  rate  given  the  range  and  signal 
strength  over  a  variety  of  system  configurations. 


II.  MINI-RANGER  III  POSITIONING  SYSTEM 


A.  SYSTEM  COMPONENTS 

A  good  understanding  of  the  components  of  the  MRS  III 
(Table  III)  and  their  operation  will  be  helpful  to  the 
reader. 


TABLE  III 


MINI-RANGER  COMPONENTS  AND  SERIAL  NUMBERS 


System  Component 
Range  Console 
RT  Unit 

Signal  Strength  Counter 
Signal  Strength  Counter 
Transponder  Code  0 
Transponder  Code  1 


Serial  Number 
(RA-4),  B0269 
(RA-5),  B1108 
No.  1  110 

No.  2  111 

C1789 

C1883 


1 .  Receiver-Transmitter  Assembly 

The  term  universal  station  refers  to  the  assembly 
that  houses  either  an  RT  unit  or  a  reference  station  (Fig. 
2.1).  A  universal  station  contains  a  motherboard,  power 
supply,  modulator,  and  the  RF  assembly.  An  RT  unit  is  a 
universal  station  with  a  video  interface  assembly  and  a 
6-db  antenna.  The  receiver- transmitter  assembly  ( RT  unit) 
consists  of  a  radar  transmitter,  a  radar  receiver,  and  a 
video  interface  assembly.  Due  to  the  complexity  of  the 
electronics,  technical  descriptions  of  the  system's  compo¬ 
nents  and  their  operation  are  taken  directly  from  the  MRS 
III  Positioning  System  Maintenance  Manual. 
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The  encoded  pulse  modulation  signal  from  the  range  con¬ 
sole  passes  through  the  video  interface  assembly  to 
the  transmitter.  The  resulting  coded  RF  signal  passes 
through  the  circulator  to  the  antenna.  At  the  same 
time,  the  transmitter  causes  a  video  signal  to  pass 
through  the  video  interface  assembly  to  the  range  con¬ 
sole  where  it  becomes  a  start  signal  for  the  range 
counters.  Meanwhile,  the  interrogation  signal  is  re¬ 
ceived,  processed,  and  a  reply  is  generated  by  the 
selected  reference  station  assembly.  The  reply  signal 
is  received  by  the  RT  unit  antenna  and  is  passed 
through  the  circulator  to  the  receiver.  The  resulting 
video  then  passes  through  the  video  interface  assembly 
and  the  interconnect  cable  to  the  decoder  assembly  in 
the  range  console.  There  it  is  decoded  and  becomes 
the  stop  signal  for  the  range  counters.  The  receiver- 
transmitter  operates  from  a  +28  and  -12  volt  DC  supply 
provided  by  the  range  console.  The  receive  and  trans¬ 
mit  frequencies  are  factory  preset  to  standard 
frequencies.  [Ref.  10:p.  7] 

2.  Shore  Reference  Stations 

An  RT  unit  is  located  on  the  survey  launch  and  in¬ 
terrogates  one  or  more  reference  stations  located  ashore. 

A  reference  station  is  a  shore-based  universal  station  that 
is  located  over  a  known  geodetic  control  station. 

A  reference  station  contains  a  16-code  decoder  as¬ 
sembly  and  a  13-db  antenna.  The  description  of  a  reference 
station's  operation  is  helpful  in  understanding  possible 
error  sources.  An  interrogation  signal  from  the  RT  unit  is 
coupled  through  the  reference  station  receiver  to  the  16- 
code  decoder  assembly.  If  the  proper  f irst-to-second  pulse 
and  second-to-third  pulse  spacing  is  detected,  the  decoder 
will  respond  with  a  two-pulse  signal  that  is  coupled  to  the 
modulator,  causing  a  reply  to  be  transmitted  back  to  the  RT 
unit.  The  f irst-to-second  pulse  spacing  is  given  in  micro¬ 


seconds  and  corresponds  to  a  given  code  value.  The 


second-to-third  pulse  spacing  is  either  50  or  54  microseconds 
A  reference  station  must  be  set  to  the  code  that  has  been  se¬ 
lected  on  the  range  console. 

3 .  Antennae 

There  are  three  different  types  of  antennae  for  use 
with  the  MRS  III  (Figs.  2.2,  2.3,  and  2.4).  The  antennae 
are  classified  by  their  purpose  which  is  related  to  the  out¬ 
put  gain  and  radiation  pattern  characteristics.  Elevation 
and  azimuth  patterns  refer  to  the  vertical  and  horizontal 
radiation  angles  of  an  antenna  (Fig.  2.5). 

The  medium-gain  antenna  (Fig.  2.2)  has  a  gain  of  13 
db  and  is  the  antenna  most  often  used  for  hydrographic  sur¬ 
vey  operations  and  base-line  calibrations.  The  horizontal 
azimuth  pattern  of  the  medium-gain  antenna  is  80  degrees  and 
its  vertical  elevation  pattern  is  15  degrees. 

The  high-gain  antenna  (Fig.  2.3)  has  a  gain  of  19  db 
and  offers  the  user  a  greater  operating  range.  The  high- 
gain  antenna  has  a  vertical  pattern  of  only  5  degrees. 

The  narrower  vertical  pattern  focuses  the  radiated  energy  to 
increase  the  operating  range  of  a  shore  station.  The  level¬ 
ing  of  a  high-gain  antenna  is  more  critical  then,  due  to  its 
narrower  vertical  pattern. 

The  omnidirectional  antenna  (Fig.  2.4)  has  a  gain  of 
6  db  and  is  used  on  the  survey  vessel  with  the  RT  unit.  It 
can  also  be  used  on  shore  reference  stations  located  on  is¬ 
lands  or  prominent  headlands  where  360-degree  azimuth 
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coverage  is  desired.  The  vertical  pattern  is  15  degrees. 


An  interesting  fact  jt  this  type  of  antenna  is  its  mis¬ 
leading  name.  "An  omni  antenna  will  receive  signals  from 
all  directions.  .  .  .  The  catch,  however,  is  that  it  does 
not  receive  equally  well  from  all  directions.  The  omni  an¬ 
tenna  has  a  preferred  direction."  [Ref.  2:p.  8] 

Previous  research  has  shown  that  the  90/270-degree 
directions  are  the  preferred  ones.  Unfortunately,  the  an¬ 
tenna  is  usually  installed  with  the  0/180-degree  axis  point¬ 
ing  toward  the  work  area.  [Ref.  10:p.  8]  The  January  1984 
calibrations  were  performed  using  the  omnidirectional  an¬ 
tenna  with  the  0/180-degree  direction  pointing  toward  the 
other  antennae. 

4 «  Range  Console 

The  range  console  (Fig.  2.6)  consists  of  a  coder 
assembly,  range  control  assembly,  range  counter  display  as¬ 
sembly,  motherboard  assembly,  and  a  power  supply  assembly. 

The  power  supply  assembly  provides  operating  voltages  for 
both  the  range  console  and  the  RT  unit.  A  description  of 
the  operation  of  the  range  console  follows.  The  coder  as¬ 
sembly  generates  the  transmitter  pulse  modulation  signal. 

When  a  start  signal  is  received  from  the  receiver-transmitter 
assembly,  the  range  counter  starts  to  count.  The  reference 
station  assembly  reply  video  is  decoded  by  the  decoder  and 
becomes  the  stop  signal  for  the  range  counter.  After  five 
consecutive  reference  station  assembly  replies,  the  range 

41 


/.V 

/  V 


y: 


V-- 


■V 

Vs/ 


?.v, 

.-  V 
/  V 


-  -S 


l 


ft 


,*  \*  *>  V  % 


counter  control  commands  the  front  panel  range  counter  dis¬ 
play  to  refresh  the  range  readout.  The  read  information,  in 
parallel  binary-coded  decimal  (BCD)  format,  also  is  available 
at  a  rear  panel  connector  for  use  by  peripheral  equipment 
such  as  printers  or  computers.  The  range  console  interro¬ 
gates  each  of  the  two  reference  station  assemblies  alternate¬ 
ly.  Two  thumbwheel  switches  on  the  front  panel  select  two 
reference  stations  by  code.  Four  distinct  reference  station 
assembly  codes  may  be  used,  two  at  a  time,  in  a  standard  MRS 
III.  An  option  to  this  allows  selection  from  any  of  16  codes 
rather  than  four,  to  improve  the  versatility  of  the  system 
should  other  users  be  present  in  the  work  area. 

B.  MRS  III  SYSTEM  OPTIONS 
1 •  Range  Average 

An  optional  feature  of  the  range  console  included  in 
all  NOS  MRS  Ill's  is  the  range  average  switch.  The  range 
average  switch  performs  internal  averaging  of  the  received 
range  rates  at  one  of  five  accuracy  levels.  This  option 
helps  the  system  produce  more  repeatable  readings. 

A  range  average  switch  is  located  on  the  front  panel 
of  the  console  that  is  labeled  in  total  turn-around  cycles 
called  PRF  cycles.  The  five  settings  are  1,  5,  20,  40,  and 
75.  "in  the  OFF  position,  the  range  console  will  take  the 


normal  5  consecutive  PRF  cycles  to  produce  one  valid 

range  reading.  In  the  75  position,  the  console  takes  75  PRF 

cycles  or  15  samples  to  automatically  provide  an  average. 


HW!" 


This  is  the  same  as  the  operator  taking  15  readings  from  the 
front  indicators  and  calculating  the  average.  If  the  range 
console  misses  one  reading  in  acquiring  the  average,  the  en¬ 
tire  group  of  readings  is  not  thrown  out,  only  that  group  of 
5  in  which  the  average  was  missed."  [Ref.  3:Section  5.2.2] 
Accuracy  claims  (by  Motorola)  for  the  MRS  III  are 
based  on  a  one-sigma  standard  deviation  or  probability  of 
68.23  percent  (Appendix  A).  The  best  results  are  obtained 
by  using  a  single  survey  vessel  and  the  75  range  average 
setting.  "If  more  than  one  range  console  is  used  in  the  sys¬ 
tem,  some  erratic  readings  will  result  if  the  RANGE  AVG 
switch  is  set  above  20."  [Ref.  3:Section  3.1.4]  Higher 
settings  require  more  update  time.  This  is  important  if 
multiple  users  need  to  share  shore  stations.  To  determine 
whether  range  averaging  is  significant  during  a  base-line 
calibration  of  a  system,  the  January  1984  calibration  data 
were  acquired  using  settings  of  20,  40,  and  75. 

2 .  Signal  Strength  Indicator 

The  signal  strength  indicator  is  mounted  atop  the 
range  console  (Fig.  2.7).  The  ranges  shown  on  the  console 
correspond  to  signal  strength  values  displayed  on  the  indi¬ 
cator.  A  brief  explanation  of  the  counter  will  aid  the 
reader  in  understanding  the  range  error/signal  strength 
relationship. 

The  video  processor  inside  the  RT  unit  converts  the 
received  pulse  amplitude  to  a  standard  amplitude  proportional 
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to  the  amplitude  of  the  received  pulse.  This  pulse  is  then 
used  to  start  and  stop  a  counter  clock  that  generates  a 
relative  strength  value  in  integer  form  [Ref.  3:Section  5.2.3] 

Typically  the  signal  strength  will  vary  from  4  to 
25  for  the  MRS  III.  NOS  counters  are  also  equipped  with  audio 
alarms  that  are  activated  when  signal  strengths  fall  below  a 
preset  value.  When  the  range  console  receives  a  signal  that 
has  undergone  attenuation,  it's  degraded  amplitude  is  con¬ 
verted  to  a  Relative  value,  displayed  on  the  indicator. 

Signal  strengths  can  be  used  with  the  alarm  option 
to  identify  weak  areas  of  coverage.  NOS  signal  strengths 
are  recorded  directly  on  the  data  tape.  The  next  step  is  to 
use  these  recorded  signal  strengths  to  extract  a  range-error 
corrector  in  near-real  time. 


III.  EXPERIMENT  DESIGN  AND  DATA  REDUCTION 


A.  OBJECTIVES 

The  main  objective  of  this  thesis  is  to  generate  an  MRS 
III  corrector  curve  that,  when  given  a  range  rate  and  corre¬ 
sponding  signal  strength  as  input,  will  yield  an  accurate 
range  corrector  that  can  be  applied  to  the  range  rate  in 
near-real  time.  The  corrector  curve  should  be  characteris¬ 
tic  of  all  possible  MRS  III  equipment  configurations  and  con 
sole  settings.  To  obtain  such  a  curve,  the  thesis  addresses 
several  related  topics: 

Range  Average  was  examined  for  its  significance  in  de¬ 
fining  the  calibration  curve. 

The  significance  of  the  three  different  types  of  antennae 
on  the  base-line  calibration  was  investigated. 

The  significance  of  base-line  distance  was  investigated 
in  hopes  of  finding  an  optimum  base-line  calibration 
length . 

An  equipment  configuration  that  is  most  characteristic  for 
an  MRS  III  unit  was  identified. 


B.  ASSUMPTIONS 

Several  assumptions  were  made  when  performing  the  field 

work  and  analyzing  the  data: 

Range  errors  were  not  related  to  power  supply  voltages 
since  the  units  furnished  by  RAINIER  were  equipped 
with  "cut-off"  switches  that  deactivate  the  transponders 
when  the  voltage  is  less  than  22  volts.  Therefore, 
errors  due  to  power  fluctuations  were  assumed  negligible. 

The  horizontal  and  vertical  pointings  of  the  reference 
station  antennae  were  considered  to  be  negligible  since 
each  code/reference  station  was  pointed  in  the  same 


manner.  The  antennae  were  pointed  visually  for  the  three 
shorter  base  lines.  For  the  two  longer  base  lines,  a 
radio  was  used  to  tell  the  reference  station  operator 
when  optimum  signal  strengths  were  observed. 

The  tellurometers  that  were  used  to  measure  the  base 
lines  were  accurate  to  plus  or  minus  0.015  meters. 

Thus,  the  lines  were  measured  to  submeter  accuracy 
so  reference  distances  introduced  no  significant 
error  into  the  data. 

Artificially  attenuating  the  incoming  signals  to  the 
RT  unit  adequately  represented  real  conditions  en¬ 
countered  in  the  field.  Though  this  assumption  is 
not  necessarily  valid,  the  artificial  attenuator 
offers  the  only  practical  means  of  simulating  the 
attenuation  caused  by  signals  traveling  long  distances. 

Internal  drift  of  the  MRS  Ill's  during  the  calibrations 
was  assumed  to  be  negligible.  Each  day  a  warm-up  time 
of  30  minutes  was  observed  before  the  first  calibration. 
Throughout  the  day,  the  units  were  not  powered  down  and 
antenna  changes  were  affected  by  momentarily  halting 
the  interrogation  process  with  the  "Hold  Display"  button 
on  the  front  of  the  range  console.  Morning  warm-up 
times  were  within  the  prescribed  time  as  recommended  in 
the  NOS  Operation  and  Installation  Manual.  Thus,  range 
errors  due  to  temperature  changes  within  the  magnetron 
of  the  MRS  III  system  were  also  considered  insignifi¬ 
cant  since  ambient  temperatures  during  the  January 
calibration  were  invariably  within  two  to  three  degrees. 
[Ref.  3] 


C.  DATA  ACQUISITION  AND  PROCESSING 

The  calibration  field  work  for  the  thesis  was  performed 
January  27  to  30,  1984,  across  Puget  Sound  and  Shilshole  Bay, 
Washington  (Fig.  3.1),  Base-line  calibrations  were  performed 
over  five  different  distances  using  the  station  at  West  Point 
Lighthouse  as  a  pivot  point  (Table  IV). 

1 •  Tellurometer  Measurements 

Prior  to  beginning  the  actual  static  calibrations, 
each  base  line  was  measured  to  submeter  accuracy  using  a 


TABLE  IV 


BASE-LINE  ENDPOINTS  AND  RESPECTIVE  LENGTHS 


Base  Line 

A 

B 

C 

D 

E 


(From  West  Point  Lighthouse) 


End  Point 


Length ( m ) 


Sewage  Gate 
Shilshole  Marina 
Meadow  Point 
Point  Monroe 
Point  Jefferson 


1061.17 

2417.48 

4083.10 

7016.77 

9861.08 


Tellurometer  model  CA-1000  electronic  distance  measuring 
instrument . 

Since  the  MRS  III  range  rates  were  slope  distances, 
no  slope  to  horizontal  distance  computations  were  necessary. 
However,  due  to  the  tellurometers  having  a  preset  value  for 
the  refractive  index  in  air  of  radio  waves  equal  to  1.000325, 
weather  data  were  taken  concurrently  to  correct  this  fixed 
value  for  variations  in  wet-  and  dry-bulb  temperatures  and 
for  barometric  pressure.  These  same  weather  data  were  input 
into  a  Hewlett-Packard  9815S  calculator.  Meteorological 
correctors  were  performed  using  an  NOS  geodetic  program 
(EDM-03/821201,  option  A). 

The  Tellurometer  observation  forms  ( NOAA  Form  76- 
61A)  for  the  five  base  lines  include  the  output  tape  gener¬ 
ated  by  the  HP  9815S  calculator  (Appendix  C  ). 

These  submeter  results  were  used  as  reference  dis¬ 
tances  for  the  five  base  lines  and  subsequently  in  the  de¬ 
termination  of  average  range  errors  to  0.1  meter. 


Data  from  the  January  1984  calibrations  were  manu¬ 
ally  logged  into  data  files  on  the  Naval  Postgraduate  School's 
IBM  3033  computer.  Ten  readings  per  signal  strength  were  re¬ 
corded  for  each  code  (codes  0  and  1).  Each  data  file  con¬ 
tained  a  header  record  that  identified  the  base  line,  code, 
range  average  setting,  antenna  type,  and  time  of  calibration. 

The  data  were  printed  in  a  readable  form  and  checked 
for  correctness.  The  data  were  also  corrected  for  a  system¬ 
atic  offset  of  16  meters  that  was  not  adjusted  or  zeroed. 
Instead  of  adjusting  the  console,  the  l’6-meter  offset  was 
applied  to  all  the  data. 

After  verifying  the  correctness  of  the  data,  a 
FORTRAN /DISSPLA  [Ref.  11]  computer  program  STAT  (Appendix  D) 
computed  mean  rates,  mean  range  errors,  and  variances  for 
both  codes. 

A  "calibration  set"  is  herein  defined  as  a  distinct 
set  of  signal  strengths  and  corresponding  range  errors  for  a 
specific  configuration  of  code  (0  or  1),  antenna  type,  range 
average  setting,  and  base  line. 

Program  SSPLOT  (Appendix  D)  plotted  mean  range  error 
versus  signal  strength  and  standard  deviation  for  each  mean 
error  (Appendix  B).  Standard  deviation  is  displayed  as  a 
vertical  bar  for  its  respective  range  error. 


Plots  of  range  error  versus  signal  strength  were 
generated  to  assess  the  significance  of  range  average  and 


antenna  type  to  the  overall  calibration  curve  for  each  code. 
Program  SCATTER  (Appendix  D)  generated  scatter  diagrams 
(Appendix  B)  of  mean  range  errors  for  each  of  the  two  codes. 

Program  WALK  (Appendix  D)  generated  plots  of  mean 
range  error  versus  base-line  distance  using  signal  strengths 
as  markers.  These  plots  (Figs.  3.2  and  3.3)  show  that  lower 
signal  strengths  consistently  separate  themselves  from  the 
others . 

4 .  Functional  Curve  Fitting 

Curve  fitting  is  a  technique  that  employs  a  least 
squares  criterion  to  determine  unique  regression  curve  of 
any  degree  up  to  a  n  minus  1,  where  n  is  the  number  of  paired 
data  points.  The  fitted  curve  is  a  functional  form  or  equa¬ 
tion  from  which  range  correctors  can  be  predicted  in  real 
time.  SDS  III  would  be  capable  of  referencing  electronic 
correctors  by  two  methods.  One  method  uses  a  table  of 
stored  data  points  in  an  array;  it  is  the  simplest  but  would 
require  some  type  of  interpolation  for  intermediate  signal 
strengths.  A  second  method  uses  the  functional  form  of  the 
calibration  curve  to  extract  range  errors  for  a  given  signal 
strength. 


The  main  advantages  of  the  functional  approach  are 
the  simplicity  of  computer  computations  and  the  accounting 
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for  intermediate  signal  strengths  in  a  nonlinear,  more  accu¬ 


rate  manner.  The  major  reason  to  limit  the  time  spent  cali¬ 
brating  electronic  control  systems  is  to  maximize  data 
acquisition  time.  NOS ' s  desire  is  to  minimize  the  calibra¬ 
tion  time  by  using  faster,  noncritical  daily  system  checks. 

If  the  calibration  procedures  were  automated,  a  greater  num¬ 
ber  of  data  could  be  sampled  and  several  base  lines  or  equip¬ 
ment  configurations  could  be  investigated.  The  time  saved 
by  letting  a  mini-computer  compute  the  statistics  and  re¬ 
gression  curves  directly  from  a  floppy  disk  containing  the 
raw  calibration  data  would  enable  the  operator  to  obtain  more 
calibration  data  in  much  less  time.  A  curve-fitting  program 
POLY  (Appendix  D),  obtained  from  Gerald's  book  on  numerical 
analysis,  was  applied  to  the  data  [Ref.  12:p.  492]. 

The  Gerald  program  computes  coefficients  for  up  to  a 
seventh-degree  polynomial  regression  curve,  though  only  fifth- 
degree  and  lower  were  attempted  on  the  January  calibration 
data.  Both  the  beta  function  (variance)  and  the  graphical 
presentation  suggest  that  a  third-degree  fit  was  the  highest 
polynomial  fit  obtainable.  Graphs  were  made  that  show  re¬ 
gression  curves  for  powers  of  one  through  three  along  with 
the  actual  mean  range  errors  that  were  used  to  produce  the 
coefficients  defining  the  regression  curves.  These  graphs 
are  discussed  in  Chapter  4.  Since  program  POLY  does  not 
print  the  regression  coefficients,  a  Hewlett-Packard  41C 
statistical  program  (SIGMAPOLYC)  was  used  to  compute  the 


coefficients  for  a  cubic  regression  curve  for  each  base  line 
for  the  two  codes.  All  data  and  programs  used  in  the  thesis 
have  been  recorded  on  9-track  magnetic  tape. 


IV.  OBSERVATIONS  AND  CONCLUSIONS 


A.  OBSERVATIONS 

1 .  Scatter  Plots 

Scatter  plots  (Appendix  B)  show  range  error  versus 
signal  strength  for  each  calibration  set.  By  inspection, 
base  lines  D  and  E  contributed  very  little  to  the  scatter 
plot  for  both  codes.  Range  error  differences  of  up  to  9 
meters  were  observed  at  the  two  lowest  signal  strengths. 

Base  line  B  exhibited  the  least  sample  deviation  and  covered 
the  greatest  range  of  signal  strengths.  Though  fewer  data 
points  were  available  for  the  two  longer  base  lines  D  and  E, 
the  general  shapes  of  these  curves  were  similar  to  the  three 
shorter  base  lines. 

Data  anomalies  occurred  for  base  line  D  code  0  (Ap¬ 
pendix  B)  between  signal  strengths  five  and  eight.  A  maxi¬ 
mum  difference  of  5  meters  was  observed  between  the 
high-gain  antenna  curves  and  the  other  two  antenna  types. 

The  agreement  between  the  three  antenna  types  was  typically 
1  to  3  meters  and  is  best  seen  on  the  mean  range  error 
plots . 

2 .  Mean  Range  Error  Plots 

Mean  range  error  plots  showed  differences  in  mean 

range  error  due  to  the  range  average  setting  were  not  signif 
icant  across  base  lines  or  antenna  types.  A  visual 
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inspection  of  the  data  showed  excellent  agreement  between 
different  range  average  settings. 

High-gain  antenna  data  showed  the  least  sample  de¬ 
viation  in  a  consistent  manner.  This  fact  and  the  great 
range  of  signal  strengths  over  which  one  can  calibrate  sug¬ 
gest  that  calibrating  with  the  high-gain  antenna  produces  a 
calibration  curve  representative  of  all  three  antenna  types. 
An  anomaly  in  the  high-gain  antenna  data  occurred  in  base 
line  A  (around  1,000  meters)  for  code  0  (Appendix  B);  The 
plot  shows  an  upward  shift  of  the  curve  to  range  errors  of 
3  to  5  meters  higher  than  the  data  from  the  other 
antennae.  The  same  behavior  was  not  observed  with  code  1. 
Inspection  of  the  mean  range  error  versus  base-line  distance 
plots  (Figs.  3.2  and  3.3)  point  to  a  MASS  of  signal  strength 
five  for  all  the  plots. 

3 .  Range  Error  Versus  Base-Line  Distance  Plots 

The  range  error  versus  base-line  distance  plots 
(Figs.  3.2  and  3.3)  show  a  dip  where  range  error  decreased 
when  signal  strength  increased  from  four  to  six  and  then  in¬ 
creased  for  signal  strengths  seven  and  above.  This  charac¬ 
teristic  dip  near  signal  strengths  four  and  five  was  noted 
for  all  equipment  configurations  and  distances.  Signal 
strength  four  was  separated  from  the  remaining  signal 
strengths  by  a  constant  positive  value  of  3  to  4 
meters  for  each  base  line.  The  other  range  errors  (with 
signal  strength  markers)  were  grouped  close  together.  The 
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characteristic  dip  mentioned  above  was  not  as  pronounced  on 
the  fitted  regression  curves--especially  the  cubic  fit. 

As  anomaly  occurs  on  the  base  line  B  range  error 
versus  base-line  distance  plot.  Base  line  B  was  partially 
obstructed  by  a  pier  200  yards  forward  of  the  West  Point 
endpoint.  The  above  mentioned  range  errors  may  have  been 
the  result  of  additional  attenuation  from  the  pier. 

4 .  Curve-Fitting  Results 

Regression  coefficients  were  computed  for  degrees 
one  through  four  (Table  V).  Variance  was  also  computed  in 
the  program  for  each  degree  of  fitted  polynomial  and  was 
used  as  a  numerical  indicator  of  a  curve's  "goodness  of  fit" 
(Table  VI).  The  regression  curves  were  plotted  (Figs.  4.1  and 

TABLE  V 

REGRESSION  COEFFICIENTS 


(y  =  a 

+  bx  + 

2 

cx  +  dx 

3  4 

+  ex 

) 

Degree  of 

Polynomial  Variance  a 

Code  0 : 

Coefficients 

bed 

e 

1 

-0.88 

0.88 

0.00 

0.00 

0.00 

2 

4.54 

-1.13 

0.06 

0.00 

0.00 

3 

17.38 

-5.66 

0.53 

-0.02 

0.00 

4 

17.34 

-5.66 

0.53 

-0.02 

0.00 

Code 

1: 

1 

0.82 

-0.02 

0.00 

0.00 

0.00 

2 

3.85 

-0.70 

0.03 

0.00 

0.00 

3 

13.34 

-4.06 

0.38 

-0.01 

0.00 

4 

13.34 

-4.06 

0.38 

-0.01 

0.00 
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TABLE  VI 


COEFFICIENTS  FOR  CUBIC  REGRESSION  CURVES 

2  3 

(y  =  a  +  bx  +  cx  +  dx  ) 


Base  Line 

VARIANCE 

a 

b 

c 

d 

Code  0 : 

A 

0.53 

18.56 

-7.47 

0.84 

-0.03 

B 

0.64 

18.14 

-5.36 

0.51 

-0.02 

C 

0.69 

22.02 

-7.86 

0.82 

-0.03 

D 

0.55 

53.44 

-21.08 

2.60 

-0.10 

E 

0.73 

34.55 

-15.82 

2.27 

-0.10 

Code  1 : 

A 

0.33 

9.11 

-2.94 

0.28 

-0.01 

B 

0.52 

15.01 

-4.22 

0.40 

-0.01 

C 

0.73 

25.49 

-8.97 

0.95 

-0.03 

D 

0.88 

66.84 

-27.09 

3.51 

-0.15 

E 

0.80 

42.67 

-18.77 

2.63 

-0.12 

4.2)  with  the  mean  range  error/signal  strength  pairs  used  as 
input  for  creating  each  higher  degree  fit.  Gerald's  book 
[Ref.  12:p.  437],  describes  the  variance  as  an  indicator  for 
deciding  what  degree  of  polynomial  should  be  used.  As  higher 
degree  polynomials  are  used,  the  deviations  of  the  points 
from  the  fitted  curve  will  decrease  until,  when  the  degree 
of  the  polynomial  n  equals  n  minus  1,  there  is  an  exact 
match  (assuming  no  duplicate  data  at  the  same  x-value). 


Variances  for  the  four  higher  polynomial  regression  curves 
indicated  that  the  cubic  fit  was  the  best.  One 
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increases  the  degree  of  approximating  polynomial  so  long  as 
there  is  a  decrease  in  the  variance.  Both  codes  had  lowest 
variances  for  the  third-degree  curve.  Knowing  this,  cubic 
regression  coefficients  were  computed  for  each  base  line  for 
the  two  codes  (Table  VII). 

The  coefficients  for  the  base  lines  were  computed  us¬ 
ing  a  Hewlett-Packard  41C  calculator  and  regression  program 
(SIGMAPOLYC) . 

The  good  agreement  between  cubic  fits  at  base  line  B 
for  each  code  and  their  respective  overall  curves  suggests 
that  a  distance  of  about  2,500  meters  is  an  optimum  range  over 
which  to  calibrate  (using  endpoints  at  sea  level).  The  cubic 
regression  curves  for  this  distance  adequately  depict  range 
error  versus  signal  strength  for  all  five  base-line  distances 
and  three  types  of  antennae. 

B.  CONCLUSIONS 

The  optimum  base-line  length  for  obtaining  a  representative 
calibration  with  both  calibration  endpoints  at  sea  level  was 
closest  to  base  line  B  (2,500  meters)  for  both  codes  0  and  1. 

At  present,  NOS  requires  the  base-line  length  to  be  between 
700  and  3,000  meters. 

The  high-gain  type  antenna  varied  slightly  more  than  the 
other  two  antenna  types.  High-gain  calibration  errors  agreed 
to  within  1  to  2  meters  with  the  medium-gain  and  omnidirec¬ 
tional  antenna  results  in  all  but  6  of  the  90  calibration  sets. 


high-gain  antennae  as  shown  by  the  code  0,  base  line  A  range- 
error  plots  (Appendix  B).  The  anomalies  suggest  that  cali¬ 
brating  with  a  high-gain  antenna  at  too  close  a  distance 
(perhaps  less  than  1,000  meters)  may  result  in  erroneous 
correctors . 

C.  RECOMMENDATIONS 

Study  the  concept  of  rejecting  data  below  MASS  and  fitting 
the  remaining  data  to  a  linear  curve. 

Investigate  spline  curve  fitting  or  the  resolving  of  a 
calibration  curve  into  linear  and  exponential  components  be¬ 
fore  deciding  upon  future  calibration  procedures. 

Give  consideration  to  mathematically  determining  the  MASS 
of  a  curve  prior  to  fitting.  Further  statistical  tests 
could  then  be  applied  to  data  below  the  MASS  to  determine 
whether  it  should  be  rejected. 

Plot  mean  range  error  versus  signal  strength  for  each 
base  line  and  for  the  cumulative  mean  range  errors  for  each 
of  the  two  codes,  thus  showing  which  base  line's  data  in¬ 
fluences  the  overall  mean  curve  the  most.  A  comparison 
should  be  made  using  a  single,  mean  corrector  MRS  III  code, 
as  is  presently  done,  with  the  predictive  results  obtained 
from  curve  fitting.  The  comparison  should  consider  the 
allowable  positioning  errors  for  different  survey  scales. 

Review  present  base-line  calibration  procedures  with 
respect  to  the  findings  of  this  thesis.  Specifically,  eval¬ 
uate  the  optimum  distance  of  2,500  meters  with  a  range 
average  setting  of  75  and  a  high-gain  antenna. 


MINI-RANGER  III  SYSTEM  SPECIFICATIONS 


Overall  System  Power  Requirements 


Range  Console 

Input  Voltage 

Power  Consumption 

El  Unii 

Input  voltage 

Power  Consumption  (at  27  VDC) 
No  interrogation  (  standby) 
Normal  interrogation 
TpapgpQpdev 

Power  Consumption  (at  27  VDC) 
No  Interrogation  ( standby) 
Normal  Interrogation 


105  to  125  volts  AC 
12  to  32  volts  DC 
60  watts  AC 
40  watts  DC 

22  to  32  volts  DC 

9  watts 
16  watts 

22  to  32  volts  DC 

8  watts 
15  watts 


Transmission  Frequencies 


fiX  Unit 


Shore 


Transmit 


Receive 


Transmit 


Receive 


5570  MHz 


5480  MHz 


5480  MHz 


5570  MHz 


Number  of  Users 

With  range  averaging  and  multi-user  options  installed. 
Motorola  recommends  no  more  than  five  range  consoles  operate 
in  a  project  area.  Although  the  consoles  may  have  different 
codes,  they  will  still  lock  each  other  out  when  transmitting 
to  a  shore  station. 


Accuracy  of  the  System 

The  following  figures  are  based  on  one-sigma  standard 
deviation.  This  means  that  68.23  percent  of  all  readings 
will  fall  within  the  quoted  accuracy. 


Range  Average  Setting 


Accuracy  (in  meters) 


3.  00 


2.  12 


1.  15 


1.  06 


0.  78 


Antennae  and  Beam  Characteristics 


;iver-Transmitter  (El  uiiifc) 

Gain  6  db 

Horizontal  360  degrees  (omnidirectional) 

Vertical  25  degrees 

£  irangRanflar 

Medium-Gain  Antenna 

Gain  13  db 

Horizontal  80  degrees 

Vertical  15  degrees 

High-Gain  Antenna 

Gain  19  db 

Horizontal  80  degrees 

Vertical  5  degrees 

Omnidirectional  Antenna 

Same  characteristics  as  the  RT  unit  antenna 


Range  of  System  Operation 


The  system  is  line-of-sight.  Though  Motorola's  opera¬ 
tion  manual  claims  a  maximum  range  of  37,000  meters,  this 
distance  is  undoubtedly  achieved  by  increasing  the  elevation 
of  either  or  both  the  RT  unit  and  shore  transponder 
antennae. 


Antenna  Type 
Medium-Gain 
High-Gain 
Omnidirectional 


Maximum  Range  (In  meters) 

32,180 

64,360 

unknown  but  less  than  Medium-Gain. 


Motorola  Options  used  by  NOS 
Option  Description  Option  Number 

Signal  Strength  2 

Range  Averaging  4 

Multiuser  5 

Byte  Test  Module  11 


30-Minute  Standby 


15 


APPENDIX  B 


RANGE  ERROR  AND  STATISTICAL  PLOTS 

The  following  graphs  are  (1)  range  error  versus  signal 
strength  (including  standard  deviations)  for  each  calibration 
set  (pages  72  through  101)  and  (2)  scatter  diagrams  of  mean 
range  error  (for  each  signal  strength)  versus  signal  strength 
for  each  of  two  MRS  III  codes  (pages  102  through  113). 
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Tv/r_  BASELINE  C 
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.  CODE  NO.  0 


SIGNAL  STRENGTH 
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...MEDIUM  CODE  NO.  0 


SIGNAL  STRENGTH 
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SIGNAL  STRENGTH 
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APPENDIX  C 


TELLUROMETER  DATA 

The  following  observation  forms  show  measurements  of  the 
five  base  lines  obtained  with  the  Tellurometers .  All  dis¬ 
tances  were  corrected  for  meteorological  effects  as  shown  on 
the  attached  HP  9815S  calculator  tapes. 


HOAA  -ORM76-6U 
(12—731 

U.S.  DEPARTMENT  OF  COMMERCE  1 
NATIONAL  OCEANIC  AND  ATMOSPHERIC  ADMINISTRATION  ' 

TELLUROMETER  OBSERVATIONS 

MODEL  CA1000 

STATE 

IVA. 

ilocality 

PU6ET  6D0KJD 

DATE 

^9  -UN  04 

STAT.ON 

(3/4SEUN£  A 

INST.  NUMBER 

1305 

HEIGHT  OF  INST. 

MO 

TO  STATION 

WEST  POIKIT 

INST.  NUMBER 

1006 

HEIGHT  OF  INST. 

170 

PROJECT  NUMBER 

CHIEF  OF  party 

OBSERVER 

BFH 

METEOROLOGICAL  DATA 

FIELD  READINGS 

T|*E  TO0  Twb 


REMOTE  K5-0H2.S 


s“"  87-0  Sfc.S 


a  \0 

3 

c 

SI 

A 

D 

A 

E 

METERS  FEET  X  0  3048006096 


,  meteorological 

i 

CORRECTION 

fixed  instrument 

CORRECTION 

CORRECTED 

slope  distance 


- o • oo& 

-  0.065 

10C?L  I7S 
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STATE 

VjA- 

LOCALITY 
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Z9  Ja W  SA 

STAT,0fSHlL5HOLE  PIER,) 
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INST.  NUMBER 

\Qo& 

HEIGHT  OF  INST. 

1.46  rv\ 

TO  STATION 

IA/E5T  POIK)T 

INST.  NUMBER 

IS05 

HEIGHT  OF  INST. 

1.77  m 

PROJECT  NUMBER 

CHIEF  OF  PARTY 

OBSERVER 

ELl-4 

meteorological  data 

FIELD  READINGS 

T  T 

OB  WB 


MASTER  M7.0  \A5.0  39.45 


SET  1 


■  ]o  |7 

c  Z5__ 
0  L£3 
e  3 


METERS  FEET  X  0.3040006096 

DISTANCE  IMefers) 

r 

i 

i 

1 

2417 .492 

METEOROLOGICA  l 

CORRECTION 

1 

1 

1 

-o,o  II 

FIXED  INSTRUMENT 

CORRECTION 

“  1- 

1 

1 

1 

'  0.004  YY\ 
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INST.  NUMBER 

1805 

HEIGHT  OF  INST. 

1  .T7 1 n 

PROJECT  NUMBER 
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&lM 

METEOROLOGICAL  DATA 

FIELD  READINGS 

1  1 

SET  1 

SET  2 

=  pk  =  -  r .  ' 

-  F  M =  -8,0;  : 

: i  ‘  !j  F  hN  •  : 

=  .  o 
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DISTANCE  IMercrs  I 


I  METEOROLOGICAL 
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—  O .  0  35 
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TIME 

T 

D8 

T 

WB 

P 

MASTER 

50.0 

49.0 

3M5 
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43.0 

41.5 

30.36 

SUM 

-D 

CD 

b 

60.21 

MEAN 

49.0 

49.2 

30.40 
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_ -SET  1 

8  IZ|3| 

c  3  O 


o_z 
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_ 1.25  m 

HEIGHT  OF  INST  . 
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APPENDIX  D 


COMPUTER  PROGRAMS 


Program 


Description 


Page 


METOBS  Uses  DISSPLA  subroutines  to  generate  plots 
of  temperature  and  barometric  pressure 
versus  time. 

RAWRITE  Reads  Mini-Ranger  rates  from  a  file; 

corrects  the  rates  for  a  16-meter  offset; 
and  writes  the  corrected  rates,  header 
information,  and  signal  strengths  to 
another  file. 

STAT  Lists  the  means  and  variances  of  both 

observed  rates  and  computed  range  errors 
of  Mini-Ranger  calibration  data. 

SSPLOT  Uses  DISSPLA  subroutines  to  generate  plots 
of  mean  range  error  (and  standard 
deviation)  versus  signal  strength. 

WALK  Uses  DISSPLA  subroutines  to  generate  plots 

of  mean  range  error  versus  base-line 
distance . 

POLYFIT  Fits  polynomial  curves  for  up  to  degree 
seven  to  the  mean  range  error  values. 

POLY  Uses  DISSPLA  subroutines  to  generate  plots 

of  range  errors  predicted  from  polynomial 
regression  curves. 

BASELINE  Uses  DISSPLA  subroutines  to  generate  plots 
of  range  error  versus  signal  strength  for 
each  base  line. 


v*v*v 

V.V/JI 

*  *  '  J>  '  i 
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PROGRAM:  METEOROLOGICAL  OBSERVATIONS  VS  TIME  (IN  HRS.) 

AUTHOR:  ERUCE  F.  HILLARD 

DATE:  AFRIL  1984 

PROGRAM  DESCRIPTION:  THIS  FORTRAN  PROGRAM  USES  DISSPLA 

SUBROUTINES  TO  GENERATE  PLOTS  OF 
TEMFERATURE  AND  BAROMETRIC  FRESSURE 
VERSUS  HOURS. 

INPUT:  INFUT  COMES  FROM  DATA  STATEMENTS. 

OUTPUT:  THE  RESULTANT  FLOT  CAN  BE  DIRECTED  TO  ONE  OF  THREE 

PLOTTERS:  VERSATEC  8222A,  VERSATEC  0900A,  OR  A 
TEKTRONIX  4631  HARDWARE  DEVICE. 

*;****************************  fc*^***********************************  A* 


REAL  TEMP, FRESS , TIME 

DIMENSION  TEMP  (50)  ,  PRESS ( 50 ) , IPAK(IOO) , TIME (50) 
**********  JAN.  27  DATA 

DATA  TIME/0400. ,0700. ,1000. ,1300. , 1600. ,1900. ,2200./ 
DATA  TEMP/47.0,48.0,49.0,50.0,51.0,50.0,50.0/ 

DATA  PRESS/30. 43 ,30.40,30.42,30.43, 30.40,30.44,30.44/ 

**********  jam. 28  DATA 

DATA  TIME/0400. ,1000. ,1300. ,1600. ,1900. ,2200./ 

DATA  TEMP/49.0,48.0,52.0,52.0,50.0,50.0/ 

DATA  PRESS/ 30.48,30.46,30.44,30.38,30.40,30.40/ 

**********  jam .  29  DATA 

DATA  TIME/0400.,0700. , 1000. ,1300.  ,1900.  ,2200./ 

DATA  TEMP/ 47. 0,46. 0,47. 0,47 .0,46.0,44.0/ 

DATA  PRESS/ 30. 42, 30. 42, 30. 42, 30. 40, 30. 35, 30. 35/ 

**********  jan. 30  DATA 


DATA  TIME/ 0400. ,1000. , 1300 . , 1600 . , 1900 . ,2200./ 
DATA  TEMF/43. 0,44. 0,44. 0,43. 0,44. 0,42.0/ 

DATA  FRESS/ 30 .32,30.30,30.26,30.22,30.18,30.17/ 

**********  INITIALIZE  PLOTTER  AND  SET  PAGE  SIZE. 

CALL  PRTPLT (72,6) 

CALL  COMPRS 

CALL  VRSTEC  (0,0,0) 

CALL  TEK618 

CALL  FACE  ( 1 1 . , 8 . 5 ) 

**********  ,jake  SS  AXIS  IN  INTEGERS. 


CALL  X INTAX 

**********  ORIENT  Y  AXIS  NUMBERS  HORIZONTALLY . 


onn  noon  n  n  non  nnnn  non  non  non  non  non 


CALL  YAXANG  (0.0) 

**********  SELECT  TOTAL  PLOT  AREA. 

CALL  AREA2D(3. ,5.) 

**********  SCALE  AND  LABEL  AXES  WITH  SELF-COUNTING  OPTION. 

CALL  YNAME (' TEMPERATURE  (DEGREES  F)$',100) 

CALL  XNAME  ('JANUARY  27,  HOURS$',100) 

CALL  BASALF( 'STAND1 ) 

**********  CREATE  THE  HEADING  AND  SET  THE  PRINT  TYPE. 

CALL  SWISSM 

CALL  HEAD IN  ('METEOROLOGICAL  DATA  AT  WEST  POINTS ',100 , 1 . 5 , 1 ) 
call  duplx 

**********  DEFINE  BOTH  AXES  AMD  FRAME  THE  SUBPLOT  AREA. 

CALL  GRAF (0200. ,200.0,7400. ,40.0,2.0,54.0) 

CALL  THKFRM( .02) 

CALL  FRAME 

**********  SET  UP  LEGEND  DATA  AND  TEXT  FOR  LEGEND  NAME. 

CALL  LINESP  (2.0) 

CALL  LINES  (' TEMPERATURES  1 , 1 PAR, 1 ) 

CALL  LINES  (' PRESSURES ', IPAK , 2 ) 

**********  FIND  SIZE  FOR  LEGEND. 

XW=XLEGND (IPAK , 2 ) 

YW=YLEGND ( IPAK , 2 ) 

CALL  MYLEGN( 'OBSERVATION  TYPE:$',100) 

**********  DEFINE  CURVE  TYPES  AND  EXECUTE  THE  CURVE  PLOT. 

CALL  LEGLIN 
CALL  DOT 

CALL  CURVE (TIME, TEMP, 7,1) 

**********  SET  UP  SECONDARY  Y-AXIS  AT  RIGHT  END. 

CALL  YGRAXS (30 .00 ,0.05,30.60,5. , 'PRESSURE  (IN.)S'  -100,8.0,0.0) 
CALL  DASH 

CALL  CURVE (TIME, PRESS, 7 ,1) 


**********  DRAW  THE  LEGEND. 

CALL  LEGEND ( IPAK ,2,4.50,1.00) 

**********  CLOSE  OUT  METAFILE  AND  RETURN  TO  LEVEL  ZERO. 

CALL  ENDPL(O) 

CALL  DOMEPL 


It] 
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************  ****************  A  **********  *************  *****  ******  *** 

PROGRAM :  RAWRITE 

DESCRIPTION:  THIS  PROGRAM  READS  MINI -RANGER  RATES  FROM  A  FILE 
AND  CORRECTS  THEM  FOR  A  16  METER  OFFSET.  THE 
PROGRAM  THEN  WRITES  THE  CORRECTED  RATES,  HEADER 
INFORMATION,  AND  SIGNAL  STRENGTHS  TO  ANOTHER  FILE. 
DATE:  MAY  20,  1984 
AUTHOR:  BRUCE  F.  HILLARD 

INPUT/OUTPUT:  BOTH  USE  A  FILEDEF  FOR  READING  AND  WRITING. 
SUBROUTINES:  NONE 

*  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  A  A  *  *  *  *  *  *  *  *  *  *  A  *  *  *  *  *  *  *  *  *  *  *  *  *  *  A  *  *  *  A  *  *  *  * 

CHARACTER  BL,A,B,C,D,E 

DATA  A  B  C  D  E  /  '  A  '  ,  1  B  '  'C1 , 'D1 , 'E' / 

INTEGER  I , J , N , SS , RAT . RATE , TIME ' RAVE , CODE , DATE , ANT , BASEL 
DIMENSION  BASEL (1500) . ANT( 1500) , CODE (1500) , TIME (1500) , RAVE (1500) , 
1  SS ( 1500 ), RAT (1500, 10) , RATE ( 1500 , 10  )  , 

1  LENGTH(1500) ,DATE(1500) 

REAL  LENGTH 

********** *************** ************************************* **** 


;o; 

;.\v> 


BEGIN  READING  THE  DATA  FROM  A  FILEDEF  INTO  ARRAYS. 

****************************************************************** 

1=0 

7  1=1+1 

READ (4, 201)  SS(I) , (RAT ( I , J) , J=1 , 10) , BASEL ( I ) , CODE ( I ) , ANT ( I ) , 

1  RAVE ( I ) , TIME ( I ) 

DO  100  J=l,10 

RATE (I , J)=RAT(I , J)-16 
100  CONTINUE 

IF (SS ( I ) . EQ . 99 )  GO  TO  50 
GO  TO  7 
50  CONTINUE 
N=I 

****************************************************************** 


AT  THIS  POINT,  THE  DATA  HAS  ALL  BEEN  STORED  IN  THE  ARRAYS  SS(I), 
RATE ( I , J ) , BASEL ( I ) , CODE ( I ) , ANT ( I ) , RAVE ( I ) ,  AND  TIME ( I ) . 

****************************************************************** 

****************************************************************** 


DETERMINE  LENGTH, BL,  AND  DATE  AND  BEGIN  WRITING  OUT  THE  DATA. 
****************************************************************** 


DO  160  1=1, M 


IF(SS ( I ) . EQ . 99 )  GO  TO 


IF ( BASEL ( I ) . EQ . 1 
IF  BASEL  I  .EQ. 2 
IF (BASEL ( I ) . EO . 3 
IF( BASEL ( I ) . EO . 4 
IF(BASEL(i) . EQ . 5 


. EQ . 1 )  GO 
. EQ . 2 )  GO 
. EO . 3 )  GO 
. EO . 4 }  GO 
. EQ . 5 1  GO 
LENGTH ( I )= 
BL=A 

DATE ( I )=30 
GO  TO  59 
LENGTH ( I  )  = 
BL=B 

DATE ( I )=29 
GO  TO  59 
LENGTH ( I  )  = 
BL=C 

DATE ( I )=27 
GO  TO  59 
LENGTH ( I  )  = 


500 
TO  51 
TO  52 
TO  53 
TO  54 
TO  55 
1061.24 


2417.48 


4083.10 


7016.77 
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BL=D 

DATE ( I )=28 
GO  TO  59 

55  LENGTH ( I) =9 861. 08 

RT  =  F 

DATE(I)=28 
GO  TO  59 

Q  ****************************************************************** 

C  CHECK  FOR  A  NEW  CALIBRATION  SET  (HEADER). 

Q  ****************************************************************** 

59  IF(I.EQ.l)  GO  TO  65 

60  IF(SS(I).LE.SS(I-1))  GO  TO  80 

65  WRITE ( 6 , 202 ) BL , LENGTH ( I ) , DATE ( I ) , TIME ( I ) , CODE ( I ) , RAVE ( I ) 

******************************** ********* ************************* 

TEST  FOR  THE  ANTENNA  TYPE  TO  BE  USED  IN  THE  HEADER. 
****************************************************************** 

IF( ANT( I ) . EQ . 1 )  WRITE ( 6 , 203 ) 

IF (ANT ( I ) .  EQ. 2 )  WRITE(6,204 
IF(ANT(I) . EQ.3)  WRITE(6,205) 

WRITE (6, 206) 

WRITE  6,207) 

*********  *:  ********  *******************************  ***************** 

BEGIN  WRITING  OUT  THE  DATA  ITSELF. 
****************************************************************** 

80  CONTINUE 

WRITE (6 , 208)  SS(I) , (RATE ( I , J) , J=1 , 10 ) 

160  CONTINUE 

****************************************************************** 
FORMAT  STATEMENTS: 

****************************************************************** 

201  FORMAT(I2,2X, 10(14, IX) ,11,11,11,12,14) 

202  FORMAT (//2X, 'BASE  LINE  1 ,A1 , 2X , 1  LENGTH  =  1 , F7 . 2 , IX , ' METERS  1 , 3X, 

1  'JANUARY  ',12,',  1984' ,2X, 14. IX, 'HOURS' ,/,18X, 'CODE:  ',I1,4X, 

1  'RANGE  AVERAGE  SETTING:  M2) 

203  FORMAT (22X, 'ANTENNA  TYPE:  MEDIUM  GAIN') 

204  FORMAT ( 2 2X, 'ANTENNA  TYPE:  OMNIDIRECTIONAL') 

205  FORMAT (22X, 'ANTENNA  TYFE  :  HIGH  GAIN') 

206  FORMAT(/ ,3X, 1  SIGNAL'  ,15X, 'OBSERVED  MINI-RANGER  RATES') 

207  FORMAT ( 2X , 'STRENGTH' ,2X, '1' ,4X, '2' ,4X, '3' ,4X, '4' ,4X, ' 5 ' ,4X, '6' , 

1  4X, '7' ,4X, '3' , 4X , '9' ,4X, ' 10 ' ) 

208  FORMAT (4X ,12, 5X ,10(14, IX)) 

500  CONTINUE 

STOP 


c 

C  PROGRAM :  STAT 

C  DESCRIPTION:  THIS  PROGRAM  LISTS  THE  MEANS  AMD  VARIANCES  OF  EOTH 

C  OBSERVED  RATES  AND  COMPUTED  RANGE  ERRORS  OF  THESIS 

C  MINI -RANGER  CALIBRATION  DATA. 

C  DATE:  MAY  20,  1984 

C  AUTHOR:  BRUCE  F.  HILLARD 

C  INPUT/OUTPUT:  BOTH  USE  A  FILEDEF  FOR  READING  AND  WRITING. 

C  SUBROUTINES :  NONE 

C 

£  *****************************  *************  •k'k'k’k-k-k'k-k-k'kk-k-k-k'k’k'klck-k-kTck-k 

c 

CHARACTER  BL,A,B,C,D,E 

DATA  A,B,C,D,E/'A' , 1 B 1 , 1 C 1 , 1 D 1 , ' E ' / 

INTEGER  I , J , N , SS , RAT , RATE .TIME , RAVE , CODE , DATE , AMT , BASEL 
DIMENSION  BASEL (1500 ) , ANT ( 1500 ) , CODE (1500 ) .TIME ( 1500 ) , RAVE ( 1500 ) , 
1  SS ( 1500 ), RAT (1500, 10) .RATE ( 1500 , 10 ) , XRATE  1500 ) , SUMH( 1500 ) , 

1  LENGTH(1500) .DATE (1500) , XR£RR(i500) , VRATE(1500) , 

1  VRERR( 1500 ) , SUMMM(1500) , RERR( 1500 , 10 ) ,RATSUM(1500) ,SUM(1500) 

REAL  LENGTH , RERR , RATSUM , XRATE , SUM , SUMM , SUMMM , XRERR , VRATE , VRERR 
C 

£  ****************************************************************** 

c 

C  BEGIN  READING  THE  DATA  FROM  A  FILEDEF  INTO  ARRAYS. 

C 

£  ****************************************************************** 
1=0 

7  1=1+1 

READ (4 , 201 )  SS(I) , (RAT( I , J) , J=1 , 10 ) , BASEL(I ) , CODE ( I ) , ANT ( I ) , 

1  RAVE(I) ,TIME(l) 

DO  100  J=l,10 

RATE ( I , J ) =RAT ( I , J ) - 1 6 
100  CONTINUE 

IF(SS ( I ) . EQ . 99)  GO  TO  50 
GO  TO  7 
50  CONTINUE 
N=I 

£  ****************************************************************** 

c 

C  AT  THIS  POINT,  THE  DATA  HAS  ALL  BEEN  STORED  IN  THE  ARRAYS  SS(I), 

C  RATE ( I , J ) , BASEL ( I ) , CODE ( I ) , ANT ( I ) , RAVE ( I ) ,  AND  TIME (I). 

C 

Q  ****************************************************************** 

c 

(2  ****************************************************************** 


s'. 


v.'. 


L' 


s 

s 


c 

C  DETERMINE  LENGTH, BL,  AND  DATE  AND  BEGIN  WRITING  OUT  THE  DATA. 

C 

Q  ****************************************************************** 

DO  160  1=1, N 

IF (SS ( I ) . EQ . 99)  GO  TO  500 
IF  BASEL(I) .EQ.l)  GO  TO  41 
IF ( BASEL ( I ) . EQ . 2 )  GO  TO  42 
IF( BASEL ( I ) . EQ . 3 )  GO  TO  43 
IF  BASEL(l) .EQ.4)  GO  TO  44 
IF(BASEL(l) .EQ.5)  GO  TO  45 

41  LENGTH ( I)=1061 .24 

BL=A 

DATE ( I )=30 
GO  TO  55 

42  LENGTH( I )=2417 . 48 

BL=B 

DATE ( I ) =29 
GO  TO  55 

43  LENGTH( I )=4083 . 10 


BL=C 

DATE ( I ) =27 
GO  TO  55 

LENGTH ( I ) =7016 . 77 


S.'  '.'.s'  s*  -s'  s'  s'  s'  s’ 


s'.s.Ls-\''s\s\'s's'‘s's-s''' 
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BL=D 

DATE ( I )=28 
GO  TO  55 

45  LENGTH (I )=9861 .08 

BL=E 

DATE ( I )=28 
GO  TO  55 

kkkkkkkkkkkkkkkkkkkkkkkkkk  ******  **************  ***********  ****  ***  kk 

BEGIN  STATISTICAL  CALCULATIONS. 

55  RATSUM( I )=0 . 0 
SUM( I )=0 . 0 

DO  56  J=l,10 

RATSUM ( I ) = RAT SUM ( I ) +FLOAT ( RATE (I , J) ) 

RERR( I , J) =RATE ( I ,  J ) -LENGTH ( I ) 

SUM( I )=SUM( I ) +RERR( I , J ) 

56  CONTINUE 

XRATE ( I )=RATSUM(I ) /10 . 0 
XRERR(I)=SUM(I)/10.0 
SUMM(I )=0 . 0 
DO  57  J=l,10 

SUMM( I )=SUMM( I )+( (RATE ( I , J ) -XRATE (I ) ) **2 ) 

57  CONTINUE 

VRATE ( I )=(SUKM( I ) ) /9 . 0 
SUMMM( I  )=0 . 0 
DO  58  J=l,10 

SUHHM(I)=SUMMM(I)+((RERR(I , J)-XRERR(I) )**2) 

58  CONTINUE 

VRERR ( I )  =  ( SUHMM ( I ) ) / 9 . 0 

59  IF(I.EQ.l)  GO  TO  65 

60  IF(SS(I) .LE.SS(I-l) )  GO  TO  80 

65  WRITE (7 , 202 )BL , LENGTH ( I ) , DATE ( I ) , TIME ( I ) , CODE ( I ) , RAVE ( I ) 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk  kkkkkkkkkkkkkkkkkkkkkkkkkkkk  ** 

TEST  FOR  THE  ANTENNA  TYPE  TO  BE  USED  IN  THE  HEADER. 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk  k  kkkkkkkkkkkkkkkkkkkkkk  *********** 


IF (ANT ( I ) . EQ . 1 )  WRITE (7 , 203 ) 
IF (ANT(I) .EQ.2)  WRITE  7,204 
IF(ANT(I) .EQ.3)  WRITE(7,205) 


WRITE (7 .206) 


Irk  irk  -kick  -A  Irk  "kit  -kick 


a 


BEGIN  WRITING  OUT  THE  DATA  ITSELF. 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


80  CONTINUE 

WRITE (7, 208)  SS ( I ) , XRATE ( I ) , VRATE ( I ) , XRERR( I ) , VRERR( I ) 
WRITE (7, 209)  RATSUM ( I ) , SUM ( I ) , SUMM ( I ) , SUNMM ( I ) 

160  CONTINUE 


FORMAT  STATEMENTS: 


201  FORMAT ( 1 2, 2X, 10 (1 4, IX), 1 1,1 1,11, 1 2, 14) 

202  FORMAT (//2X, 'BASE  LINE  ', A1 , 2X , 1  LENGTH  =  '  , F7 . 2 , IX , ’ METERS  1  , 3X , 
1  'JANUARY  ',12,',  1984' ,2X, 14, IX, 'HOURS' ,/,18X, 'CODE:  ’,I1,4X, 

1  'RANGE  AVERAGE  SETTING:  M2) 

203  FORMAT (22X, 'ANTENNA  TYPE:  MEDIUM  GAIN') 

204  FORMAT  22X , 1  ANTENNA  TYPE:  OMNIDIRECTIONAL') 

205  FORMAT  22X, 'ANTENNA  TYPE:  HIGH  GAIN') 

206  FORMAT (/ ,4X, 'SIGNAL' ,5X, 'AVERAGE1 ,7X, 'RATE' ,8X, 'AVERAGE' ,5X, 

1  ' RANGE  ERROR '  , / , 3X , ' STRENGTH  1  , 6X , ' RATE '  , 6X ,  ' VARIANCE  1  , 5X , 

1  'RANGE  ERROR'  ,3X, 'VARIANCE  1  ) 
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PROGRAM:  SSPLOT 
AUTHOR:  BRUCE  F.  HILLARD 

DATE:  JUNE  6,  1934 

PROGRAM  DESCRIPTION:  THIS  FORTRAN  PROGRAM  USES  DISSPLA 

SUBROUTINES  TO  GENERATE  PLOTS  OF 
MINI -RANGER  SIGNAL  STRENGTH  VERSUS 
MEAN  RANGE  ERROR  WITH  STANDARD 
DEVIATIONS. 

INPUT:  INPUT  COMES  FROM  A  DATA  FILE  REFERENCED  WITH  A 

FILEDEF. 

OUTPUT:  THE  RESULTANT  PLOT  CAN  BE  DIRECTED  TO  ONE  OF  THREE 

PLOTTERS:  VERSATEC  8222A,  VERSATEC  0900A,  OR  A 
TEKTRONIX  4631  HARDWARE  DEVICE. 


INTEGER  I, J,K,L,M,SS, TIME, RAVE, CODE, ANT, BASEL, 

1  SSA, SSI , SS2 , SS3 , P , Q , R , S , AMTA , IPAK 

DIMENSION  BASEL ( 1 500 ), ANT (1500) ,CODE(1500) , TIME (1500) , 

1  SS ( 1500 ) , XRERR (1500) , XRATE (1500) , VRERR( 1 500 ) , VRATE (1500) , 

1  XRERR1 ( 1500 ) , XRERR2 (1500) , XRERR3 ( 1 500 ), XRERRA ( 1500 )  , 

1  SSA( 1500) , SSI (1500) ,SS2( 1500) ,SS3( 1500) , RAVE (1500) ,STD1A(1500) , 

1  STD1B(1500) ,STD2A( 1500) ,STD2B( 1500) , STD3A( 1500 ), STD3B ( 1500 ) , 

1  VRERRA(1500) ,VRERR1( 1500) ,VRERR2( 1500) , VRERR3( 1500) ,ANTA(1500) , 

1  TEMPI (1500) , TEHP2 ( 1500 ) ,TEMP3(1500) ,IPAK(100) 

REAL  XRERR , XRERRA , XRERR1 , XRERR2 , XRERR3 , XRATE , VRATE , VRERR , 

1  STD1A , STD1B , STD2A , STD2B , STD3A , STD3B , VRERRA , VRERR1 , VRERR2 , VRERR3 , 
1  TEMPI ,TEHP2,TEMP3 

**********  READ  IN  THE  STAT  DATA  FROM  A  FILEDEF. 


J=0 

K=0 

L=0 

M=0 

1=0 

P=1 

0=1 

R=1 

S=1 

7  1=1+1 

READ(4,208)  SS ( I ) , XRATE ( I ) , VRATE ( I ) , XRERR( I ) , VRERR( I ) ,BASEL(I) , 
1CODE ( I ) , ANT ( I ) , RAVE ( I ) , TIME ( I ) 

IF(SS ( I ) .EQ.99)  GO  TO  50 
GO  TO  7 

50  CONTINUE 
N=I 

**********  AT  THIS  POINT,  ALL  THE  DATA  HAS  BEEN  STORED  IN  ARRAYS. 

**********  DETERMINE,  BY  COMPARISONS,  THE  DATA  TO  BE  PLOTTED. 

DO  160  1=1, N 

55  I F ( BASEL ( I ) . EQ . 5 )  GO  TO  60 

GO  TO  160 

60  IF ( CODE ( I ) . EQ . 0 )  GO  TO  70 

GO  TO  160 

70  IF ( RAVE ( I ) . EQ . 20 )GO  TO  80 

GO  TO  160 

30  SSA ( P )=SS  ( I ) 

XRERRA ( P ) =XRERR ( I ) 

VRERRA ( P ) =VRERR ( I ) 

AMTA ( P ) =ANT ( I ) 

P=P+1 


160  CONTINUE 


non  non  onn  nnn  nnn  nnonoon  n  oooonoo 


J=P-1 

c**************  at  THIS  POINT,  ONLY  THE  "BASE  LINE  A/CODE  0/RAVE=20" 
VALUES  HAVE  BEEN  ACCEPTED  INTO  NEW  ARRAYS 
SUBSCRIPTED  BY  THE  INTEGER  P. 

*************  N0W  L0AD  tHESE  VALUES  BY  ANTENNA-TYPE  INTO  THREE 
SEPARATE  ARRAYS  AND  CALCULATE  THEIR  STANDARD 
DEVIATIONS  FROM  THEIR  MEANS. 

DO  180  P=1,J 

IF(AMTA(P) .EQ.l)  GO  TO  191 
IF  ANTA  P  .EQ.2  GO  TO  192 
IF(ANTA(P) .EQ.3)  GO  TO  193 

191  SS1(Q)=SSA(P) 

KRERR1 (Q)=XRERRA(P ) 

VRERR1 (Q)=VRERRA(P ) 

STD1A(Q)=XRERR1(Q)+SQRT(VRERR1(Q) ) 

STD 1 8 ( Q ) =XRERR1 ( Q ) - SQRT ( VRERR1 ( Q ) ) 

Q=Q+1 
GO  TO  180 

192  SS2(R)=SSA(P) 

XRERR2 ( R ) =XRERRA ( P ) 

VRERR2 (R)=VRERRA(P) 

STD2A(R)=XRERR2 (R)+SQRT (VRERR2 (R) ) 

STD2B (R)=XRERR2 ( R) -SQRT (VRERR2 (R) ) 

R=R+1 
GO  TO  180 

193  SS3(S)=SSA(P) 

XRERR3 ( S ) =XRERRA ( P ) 

VRERR3 ( S ) =VR£RRA ( P ) 

STD3A(S )=XRERR3 (S )+SQRT (VRERR3 (S ) ) 

STD3B ( S ) =XRERR3 ( S ) - SQRT ( VRERR3 ( S ) ) 

S=S+1 
GO  TO  180 

180  CONTINUE 
K=Q-1 
L=R-1 
M=S-1 


*************  PERFORM  THE  NECESSARY  PLOT/SETUP  ROUTINES. 

**********  INITIALIZE  PLOTTER  AND  SET  PAGE  SIZE. 

CALL  PRTPLT (72,6) 

CALL  COtlPRS 

CALL  VRSTEC  (0,0,0) 

CALL  TEK618 

CALL  PAGE  (11. ,8.5) 

CALL  BLOWUP (1.25) 

**********  hake  ss  AXIS  IN  INTEGERS 
CALL  XINTAX 

**********  ORIENT  Y  AXIS  NUMBERS  HORIZONTALLY. 

CALL  YAXANG  (0.0) 

**********  SELECT  TOTAL  PLOT  AREA. 

CALL  AREA2D ( 8 . ,6.  ) 

**********  SCALE  AND  LABEL  AXES  WITH  SELF-COUNTING  OPTION. 

CALL  YNAMEC RANGE  ERROR  (M)  S', 100) 

CALL  XNAHE  ('SIGNAL  STRENGTHS  1  , 100 ) 

CALL  BASALF( ’STAND' ) 

**********  CREATE  THE  HEADING  AND  SET  THE  PRINT  TYPE. 


nnn  onn  nnn  nnn  nnn  nnnn  n  n  n  n  nnn  non 


CALL  SWISSM 

CALL  HEADIN  ('MINI-RANGER  III  BASELINE  CALIBRATIONS  1 , 100 , 1 . 5 , 1 ) 
**********  DEFINE  BOTH  AXES  AND  FRAME  THE  SUBPLOT  AREA. 

CALL  DUPLX 

CALL  GRAF( 18.0, 1.0,3.0,-10.0,2.0,10.0) 

CALL  THKFRM( .02) 

CALL  FRAME 

**********  SET  UP  LEGEND  DATA  AND  TEXT  FOR  LEGEND  NAME. 

CALL  LINESP  (2.0) 

CALL  LINES  ( 1  MEDIUMS IPAK, 1 ) 

CALL  LIMES  ( 1 OMNIS ' , IPAK, 2) 

CALL  LINES  ( 1  HIGHS IPAK, 3 ) 


**********  FIND  SIZE  FOR  LEGEND. 


XW=XLEGND ( IPAK , 3 ) 

YW=YLEGND (IPAK , 3 ) 

CALL  MYLEGN ( ' ANTENNA  TYPE : $ ' , 100 ) 

CALL  LEGLIN 

**********  DEFINE  CURVE  TYPES  AND  EXECUTE  THE  CURVES. 


**********  PLOT  THE  MEDIUM  GAIN  ANTENNA  VALUES. 


CALL  MARKER( 15 ) 
CALL  DOT 


CALL  THKCRV ( .025) 

DO  195  0=1, K 

TEMPI (Q)=FLOAT(SSl (Q) ) 
195  CONTINUE 


CALL  CURVE (TEMPI , XRERR1 , K , 


1) 


***********  PLOT  THE  OMNIDIRECTIONAL  ANTENNA  VALUES. 


CALL  MARKER (16) 

CALL  DASH 

CALL  THKCRV (.025) 

DO  196  R=1,L 

TEMP2 (R)=FLOAT(SS2(R) ) 

196  CONTINUE 

CALL  CURVE ( TEMP 2 , XRERR2 , L , 1 ) 


************  PLOT  THE  HIGH  GAIN  ANTENNA  VALUES. 


CALL  MARKER (17) 

CALL  RESET ( 'DASH' ) 

CALL  THKCRV ( . 025 ) 

DO  197  S=1 , M 

TEMP3(S)=FL0AT(S33(S) ) 

197  CONTINUE 

CALL  CURVE ( TEMP  3 , XRE  RR  3 , M , 1 ) 


**********  OTHER  PLOT/LEGEND  INFORMATION: 


CALL  MESSAG( 'BASELINE  ES ' , 100 , 3 . 50 , 5 . 50 ) 

CALL  MESSAG( 1  RANGE  AVERAGE:  20S 1  , 100 , 3 . 50 , 5 . 25 ) 

CALL  MESSAG( 'RANGE  AVERAGE:  40S ' . 100 , 3 . 50 , 5 . 25 ) 

CALL  MESSAG( 'RANGE  AVERAGE:  75S ' , 100 , 3 . 50 , 5 . 25 ) 

CALL  MESSAGE  CODE  NO.  OS  1 , 100 , 3 . 5 , 5 . 0 ) 

CALL  MESSAGl ' CODE  NO.  IS ' , 100 , 3 . 5 , 5 .0) 

CALL  MESSAG( 'NOTE:  VERTICAL  BARS  ARE  STANDARD  DEVIATIONS . S' , 
1  100,1.25,0.50) 


**********  DRAW  THE  legend. 


CALL  LEGEND ( IPAK ,3,1.0,3.75) 
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c  *******  PLOT  vertical  bars  (standard  DEVIATIONS), 

c 

CALL  BARWID (0 . 025 ) 

CALL  8ARPAT(16) 

CALL  VBARS (TEMPI , STD1B , STD1A ,  K) 

CALL  VBARS (TEMP 2 , STD2B , STD2A ,  L) 

CALL  VBARS (TEMP 3 , STD3B , STD3A ,M) 

C 

C  **********  READ  format  statement. 

c 

208  FORMAT (6X, 12 , 7X, F6 . 1 , 8X, F3 . 1 , 10X, F4. 1 , 11X, F3 . 1 , 3X , II , II , II , 12 , 14) 
C 

C  **********  CLOSE  OUT  METAFILE  AND  RETURN  TO  LEVEL  ZERO. 

C 

CALL  ENDPL(O) 

CALL  DONEPL 
C 
C 

STOP 

END 
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PROGRAM:  WALK 

AUTHOR:  BRUCE  F.  HILLARD 

DATE:  AUGUST  24,  1984 

PROGRAM  DESCRIPTION:  THIS  FORTRAN  PROGRAM  USES  DISSPLA 

SUBROUTINES  TO  GENERATE  PLOTS  OF 
BASE  LINE  DISTANCE  VERSUS 
MEAN  RANGE  ERROR.  SIGNAL  STRENGTH 
VALUES  ARE  USED  AS  MARKERS. 

INPUT:  INPUT  COMES  FROM  A  DATA  FILE  REFERENCED  WITH  A 

FILEDEF. 

OUTPUT:  THE  RESULTANT  PLOT  CAN  BE  DIRECTED  TO  ONE  OF  THREE 

PLOTTERS:  VERSATEC  8,222k,  VERSATEC  0900A,  OR  A 
TEKTRONIX  4631  HARDWARE  DEVICE. 

********************************************************************* 


INTEGER  I , J , K , TIME , RAVE , CODE , ANT , BASEL , SSA , SS , BASELA , B 
IN , NN , NNN , NNNN , M , MM , MMM , MMMM ,P,PP,PPP,PPPP,R,RR, RRR 
DIMENSION  BASEL ( 1500) , ANT (1500) ,CODE(1500) ,TIME(1500)  , 

1  SS(1500) ,XRERR( 1500) , XRATE ( 1500) ,VRERR( 1500) ,VRATE( 1500) , 

1  XRERRA( 1500) ,SSA(1S00) , RAVE (1500) , BASELA (1500) 

REAL  XRERR ,  XRERRA  ,  XRATE  ,  VRATE  ,  VREP.R  ,  L 
1  X4 , XX4 , X5 , XX5 , X6 , XX6 , X 7 , XX7 , X8 , XX3 , X9 , XX9 , X10 , XX10 , 

1  X11,XX11,X12,  XXI 2  ,  X13 ,  XXI 3  ,  X14 ,  XX14 ,  XI 5  ,XX15 ,  X16 ,  XX16  ,X17  ,  XXI 7  , 
1  X18 , XX18 

**********  INITIALIZE  PLOTTER  AND  SET  PAGE  SIZE. 

CALL  PRTPLT (72,6) 

CALL  COMPRS 

CALL  VRSTEC  (0,0,0) 

CALL  TEK613 

CALL  PAGE  (11. ,8.5) 


**********  ORIENT  Y  AXIS  NUMBERS  HORIZONTALLY. 

CALL  YAXANG  (0.0) 

**********  SELECT  TOTAL  PLOT  AREA. 

CALL  AREA2D ( 8 . ,6. ) 

**********  SCALE  AND  LABEL  AXES  WITH  SELF-COUNTING  OPTION. 

CALL  YNAME ( 1  RANGE  ERROR  (M)  S'  . 100 > 

CALL  XNAME  (‘BASE  LINE  DISTANCE  (M)$',100) 

CALL  BASALF( 'STAND' ) 

**********  CREATE  THE  HEADING  AND  SET  THE  PRINT  TYPE. 

CALL  HEADIN  ('BASE  LINE  DISTANCE  VS.  RANGE  ERRORS  1 , 100 , 1 . 5 , 1 ) 
**********  DEFINE  BOTH  AXES  AND  FRAME  THE  SUBPLOT  AREA. 

CALL  DUPLX 

CALL  GRAF (0.0,1000.0,11000.0,-3.0,1.0,6.0) 

CALL  THKFRM( . 02 ) 

CALL  FRAME 


**********  OTHER  PLOT/LEGEND  INFORMATION: 

CALL  MESSAG ( 1  CODE  MO.  OS '  , 100 , 3 . 0 . 5 . 25 ) 

CALL  MESSAG( 'NUMBERS  ARE  SIGNAL  STRENGTHS  S', 
1  100,2.50,5.50) 
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c 

c  **********  read  in  the  stat  data  from  a  filedef. 
c 

J=0 

1=0 

7  1=1+1 

READ( 4, 300)  SS ( I ) , XRATE ( I ) , VRATE ( I ) , XRERR( I ) , VRERR( I ) , BASEL ( I ) , 

1 CODE ( I ) . ANT ( I ) . RAVE ( I) , TIME ( I ) 

IF(SS(I) .EQ.99)  GO  TO  50 
GO  TO  7 
50  CONTINUE 
N=I 
C 
C 

c  **********  READ  FORMAT  STATEMENT. 

C 

300  FORMAT (6X, 12 ,7X, F6 .1,8X,F3.1,10X,F4.1,11X,F3.1,3X,I1,I1,I1,I2,I4) 
C 

C  **********  AT  THIS  POINT,  ALL  THE  DATA  HAS  BEEN  STORED  IN  ARRAYS. 

C 

C  **********  SORT  OUT  THE  DESIRED  CODE. 

C 

C 

DO  400  1=1, N 

IF (CODE ( I ) .NE . 0)  GO  TO  400 
J=J+1 

SSA ( J)=SS ( I ) 

XRERRA ( J ) =XRERR ( I ) 

BASELA ( J ) =BASEL ( I ) 

400  CONTINUE 
C 

c**************  K  IS  THE  TOTAL  NUMBER  OF  POINTS  TO  BE  PLOTTED. 

C 

K=J 

B=0 

C 

c**************  D0  LOOP  THROUGH  THE  DATA  ONCE  PER  BASE  LINE. 

C 

500  CONTINUE 
B=B+1 

IF (B . EQ . 6 )  GO  TO  2000 
C 

c**************  seT  all  COUNTERS  TO  ZERO. 

C 

N=Q 

X4=0 . 0 
XX4=0 . 0 
NM=0 
X5=0 . 0 
XX5=0 . 0 
NMN=0 
X6=0 . 0 
XX6=0 . 0 
NNNN=0 
X7=0 . 0 
XX7=0 . 0 
M=0 
X8=0 . 0 
XX8=0 . 0 
MM=0 
X9=0 . 0 
XX9=0 . 0 
MMM=0 
X10=0.0 
XX10=0.0 
MMMM=0 
X11=0.0 
XXI 1=0.0 
P=0 
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X12=0 . 0 
XX12=0 . 0 
PP=0 
X13=0 . 0 
XX13=0 . 0 
PPP=0 
X14=0 . 0 
XX14=0 . 0 
PPFP=0 
X15=0 . 0 
XXI 5=0.0 
R=C 

X16=0 . 0 
XXI 6=0.0 
RR=0 
XI  7=0 . 0 
XX17=0 . 0 
RRR=0 
X13=0 . 0 
XXI 8=0.0 
C 

c**************  SORT  THE  DATA  BY  BASE  LINES  AND  SIGNAL  STRENGTHS 
C 

DO  1000  J=1,K 
C 

c**************  S0RT  F0R  0NE  BASE  LINE  AT  A  TIME. 

C 

IF(BASELA( J) .ME.B)  GO  TO  1000 
C 

c **************  S0RT  the  SIGNAL  STRENGTHS. 

C 


100  CONTINUE 

IF (SSA( J ) .NE .4)  GO  TO  101 
X4=X4+XRERRA ( J ) 

M=N+1 

GO  TO  1000 

101  CONTINUE 

IF(SSA(J) .NE.5)  GO  TO  102 
X5=X5+XRERRA( J) 

NN=NN+ 1 
GO  TO  1000 

102  CONTINUE 

IF ( SSA( J ) .ME . 6 )  GO  TO  103 

X6=X6+XRERRA(J) 

NNM=NNN+1 
GO  TO  1000 

103  CONTINUE 

IF ( SSA ( J ) . NE . 7 )  GO  TO  104 
X7=X7+XRERRA( J ) 

NNNN=NNNN+1 
GO  TO  1000 

104  CONTINUE 

IF (SSA( J ) . NE . 8 )  GO  TO  105 
X8=X8+XRERRA ( J ) 

11=11+1 

GO  TO  1000 

105  CONTINUE 

IF ( SSA( J ) .ME. 9)  GO  TO  106 
X9=X9+XRERRA( J) 

MM=IIM+1 
GO  TO  1000 

106  CONTINUE 

IF (SSA( J ) . ME . 10 )  GO  TO  107 
X10=X10+XRERRA( J) 

NNN=MHM+1 
GO  TO  1000 

107  CONTINUE 

I F ( SSA ( J ) .ME. 11)  GO  TO  108 
X11=X11+XRERRA( J ) 
MI1HM=MMMH+1 


GO  TO  102 


GO  TO  103 


GO  TO  104 


GO  TO  105 


GO  TO  106 


GO  TO  107 


on 


GO  TO  1000 

108  CONTINUE 

IF (SSA( J ) .ME. 12)  GO  TO  109 
X12=X12+XRERRA( J ) 

P=P+1 

GO  TO  1000 

109  CONTINUE 

IF (SSA( J ) . NE . 13 )  GO  TO  110 
X13=X13+XRERRA( J ) 

PP=PP+1 
GO  TO  1000 

110  CONTINUE 

IF(SSA(J) .NE.14)  GO  TO  111 
X14=X14+XRERRA( J) 

PPP=PPP+1 
GO  TO  1000 

111  CONTINUE 

IF(SSA(J) .NE.15)  GO  TO  112 
X15=X15+XRERRA( J) 

PPPP=PPPP+1 
GO  TO  1000 

112  CONTINUE 

IF(SSA(J) .ME. 16)  GO  TO  113 
XX16=X16+XRERRA( J ) 

R=R+1 

GO  TO  1000 

113  CONTINUE 

IF(SSA( J) .ME. 17)  GO  TO  114 
X17=X17+XRERRA(J) 

RR=RR+1 
GO  TO  1000 

114  CONTINUE 

IF(SSA(J) .ME. 18)  GO  TO  1000 
X18=X18+XRERRA( J ) 

RRR=RRR+1 
1000  CONTINUE 
C 

c**************  COMPUTE  THE  MEAN  RANGE  ERROR  FOR  EACH  SIGNAL  STRENGTH 
AND  ASSIGN  THE  BASE  LINE  LENGTH  TO  BE  PLOTTED. 

IF(B.NE.l)  GO  TO  90 
L=1061 .24 
GO  TO  600 

90  IF(B .NE . 2 )  GO  TO  91 
L=2417 . 43 
GO  TO  600 

91  I F ( B . NE . 3 )  GO  TO  92 
L=4083 . 10 
GO  TO  600 

92  IF ( B .NE . 4 )  GO  TO  93 
L=7016 . 77 
GO  TO  600 

93  CONTINUE 
L=9861 . 08 

600  CONTINUE 
c**************  COMPUTE  THE  MEAN  RANGE  ERRORS. 

C 

XX4=X4/ FLOAT (N) 

XX5=X5/ FLOAT (MN) 

XX6=X6/ FLOAT (NNN) 

XX7 =X7/ FLOAT  NNNN) 

XX8=X8/ FLOAT (M) 

XX9=X9/ FLOAT (MM) 

XX10=X10/ FLOAT (MMM) 

XX11=X1 1/ FLOAT (MMMM) 

XX12=X12/FLOAT ( P ) 

XX13=X13/FLOAT(PP) 

XX14=X 14/ FLOAT ( PPP ) 

XXI 5=X1 5/ FLOAT ( PPPP ) 

XX16=X16/ FLOAT  R) 


XXI 7=X1 7 / FLOAT ( RR ) 

XX18=X18/FL0AT(RRR) 

C 

c**************  SET  THE  CHARACTER  HEIGHT. 

C 

CALL  HEIGHT (0.15) 

C 

c**************  PL0T  THE  VALUES  FOR  ONE  BASE  LINE. 

C 

IF(XX4 . EQ .0 . 0)  GO  TO  201 
CALL  RLINT(4,L,XX4) 

201  CONTINUE 

IF (XX5 . Ep . 0 . 0)  GO  TO  202 
CALL  RLINT (5 , L,XX5) 

202  CONTINUE 

IF (XX6 . EQ . 0 . 0 )  GO  TO  203 
CALL  RLINT (6 , L ,XX6) 

203  CONTINUE 

IF (XX7 . EQ . 0 . 0 )  GO  TO  204 
CALL  RLINT(7 , L,XX7) 

204  CONTINUE 

IF (XX8 . EQ . 0 . 0 )  GO  TO  205 
CALL  RLINT(8,L,XX8) 

205  CONTINUE 

IF(XX9 . EQ . 0 . 0 )  GO  TO  206 
CALL  RLINT (9 , L ,XX9 ) 

206  CONTINUE 

IF(XXIO.EQ.O.O)  GO  TO  207 
CALL  RLINT ( 10 ,L,XX10) 

207  CONTINUE 

IF (XXII . EQ . 0 . 0 )  GO  TO  208 
CALL  RLINT (11 , L,XX11 ) 

208  CONTINUE 

IF (XX12 .EQ . 0 .0)  GO  TO  209 
CALL  RLINT ( 12, L.XX12) 

209  CONTINUE 

IF (XX13 . EQ . 0 . 0 )  GO  TO  210 
CALL  RLINT ( 13, L.XX13) 

210  CONTINUE 

IF (XX14 .EQ . 0 . 0 )  GO  TO  211 
CALL  RLINT ( 14, L,XX14) 

211  CONTINUE 

IF (XXI 5 . EQ . 0 . 0)  GO  TO  212 
CALL  RLINT ( 1 5 , L , XXI 5 ) 

212  CONTINUE 

IF (XX16 . EQ . 0 . 0 )  GO  TO  213 
CALL  RLINT ( 16 , L , XX16 ) 

213  CONTINUE 

IF (XX17 . EQ . 0 . 0 )  GO  TO  214 
CALL  RLINT ( 17 ,L,XX17) 

214  CONTINUE 

IF ( XX18 .EQ . 0 . 0 )  GO  TO  215 
CALL  RLINT ( 18, L.XX18) 

215  CONTINUE 
C 

GO  TO  500 
2000  CONTINUE 
C 

c*************  CLOSE  OUT  METAFILE  AND  RETURN  TO  LEVEL  ZERO. 
C 

CALL  ENDPL(O) 

CALL  DONEPL 

STOP 

END 
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IMPLICIT  REAL*4  (A-H,0-Z) 

DIMENSION  X(iOO) ,Y(100) ,  C(1C0) , A( 10 , 1 1 ) .  XN(100) 

DATA  X(l) ,X(2) ,X(3) ,X(4) ,X(5) ,X(6) ,X(7) ,X(8) ,X(9) ,X(10) ,X(11) , 
X(12) , X ( 1 3 ) 


A****************************************************************** 

PROGRAM:  POLYFIT 
AUTHOR:  BRUCE  F.  HILLARD 

DATE:  SEPTEMBER  1984 

PROGRAM  DESCRIPTION:  THIS  FORTRAN  PROGRAM,  TAKEN  FROM 

GERALD'S  NUMERICAL  ANALYSIS  BOOK, 

FITS  POLYNOMIAL  CURVES  FOR  UP  TO  DEGREE 
7  TO  THE  MEAN  RANGE  ERROR  VALUES 
FOR  MINI-RANGER  CODES. 

INPUT:  MEAN  RANGE  ERRORS  ARE  INCLUDED  IN  A  DATA  STATEMENT. 

OUTPUT:  THE  RESULTS  ARE  WRITTEN  TO  A  PRINTER  OR  DISK  FILE. 


1  X(  14 i/17. 0’,  i6. 0,15. 0,14. 0,13. 0,12. 0,11. 0,10. 0,9. 0,8. 0,7.0, 

1  6. 0,5. 0,4.0/ 

DATA  Y(l) ,Y(2) ,Y(3) ,Y(4) ,Y(5) ,Y(6) ,Y(7) ,Y(3) ,Y(9) ,Y(10) , Y ( 1 1 ) , 

1  Y ( 12 ) , Y( 13 ) , 

1  Y( 14); 0.4, 1.4, 0.8, -0.3, 0.5, 0.6, 0.8, 0.5, 0.4, 0.0, -0.5, -0.9, 0.5, 

1  3.6/ 

THIS  PROGRAM  IS  USED  TO  IN  FITTING  A  POLYNOMIAL  TO  A  SET  OF  DATA. 

THE  PROGRAM  READS  IN  M  PAIRS  OF  X  AND  Y  VALUES  AND  COMPUTES 
THE  COEFFICIENTS  OF  THE  NORMAL  EQUATIONS  FOR  THE  LEAST  SQUARES 
METHOD . 

PARAMETERS  ARE: 

X , Y  ARRAY  OF  X  AND  Y  VALUES. 

N  NUMBER  OF  DATA  PAIRS 

MS ,MF  THE  RANGE  OF  DEGREE  OF  POLYNOMIALS  TO  BE  COMPUTED. 

THE  MAXIMUM  DEGREE  IS  9. 

A  AUGMENTED  ARRAY  OF  THE  COEFFICIENTS  OF  THE  NORMAL 

EQUATIONS. 

C  ARRAY  CF  COEFFICIENTS  OF  THE  LEAST  SQUARES  POLYNOMIALS. 


N=14 

MS=1 

MF=7 

C  COMPUTE  MATRIX  OF  COEFFICIENTS  AND  R.H.S.  FOR  MF  DEGREE. 

C  HOWEVER,  FIRST  CHECK  TO  SEE  IF  MAX  DEGREE  REQUESTED  IS  TOO  LARGE. 

C  IT  CANNOT  EXCEED  N-l.  IF  IT  DOES,  REDUCE  TO  EQUAL 
C  N-l  AND  PRINT  MESSAGE. 

IF(MF.LE. (N-l) )  GO  TO  5 
MF-N-1 

WRITE (6, 200)  MF 

200  FORMAT (/, 'DEGREE  OF  POLYNOMIAL  CANNOT  EXCEED  N-l1,//, 

1  'REQUESTED  MAXIMUM  DEGREE  OF  POLYNOMIAL  TOO  LARGE-REDUCED  TO1 
1  13) 

5  MFP1=MF+1 

MFP2=MF+2 

C  PUT  ONES  INTO  A  NEW  ARRAY.  THIS  WILL  HOLD  POWERS  OF  THE  X  VALUES 
C  AS  WE  PROCEED. 

DO  10  1=1, N 
10  XN ( I ) = 1 

C  COMPUTE  FIRST  COLUMN  AND  N+l  ST  COLUMN  OF  A. 

C  ROWS,  J  SUMS  OVER  THE  N  VALUES. 

DO  30  1=1 , MFP1 
A ( I , 1 )=0 . 0 
A( I ,MFP2)=0 .0 
DO  20  J=1 , N 
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A(  1 , 1 )  =  A(I , 1 )  +  XN(  J) 

A( I ,MFP2 )=  A(I ,MFP2)  +  Y(J)*XN(J) 

20  XN ( J ) =  XN(J)  *  X(J) 

30  CONTINUE 

C  COMPUTE  THE  LAST  ROW  OF  A.  I  MOVES  ACROSS  THE  COLUMNS, 

C  J  SUMS  OVER  N  VALUES. 

C 

DO  50  1=2 ,MFP1 
A(MFP1 , I )  =  0.0 
DO  40  J  =  1,N 
40  XN ( J )  =  XN( J)  *  X(J) 

50  CONTINUE 

C  NOW  FILL  IN  THE  REST  OF  THE  A  MATRIX.  I  MOVES  DOWN  THE 
C  ROWS,  J  MOVES  ACROSS  THE  COLUMNS. 

DO  70  J  =  2 , MFP1 
DO  60  I  =  1 ,MF 
60  A( I , J)  =  A(I+1 , J-l) 

70  CONTINUE 

C  WRITE  OUT  THE  MATRIX  OF  NORMAL  EQUATIONS. 

WRITE (6 , 201 )  ( (A( I , J) ,  J  =  1,MFP2),I  =  1,MFP1) 

201  FORMAT (/,9F13. 5) 

C  NOW  CALL  A  SUBROUTINE  TO  SOLVE  THE  SYSTEM.  DO  THIS  FOR  EACH  DEGREE 

C  FROM  MS  TO  MS. 

C  GET  THE  LU  DECOMPOSITION  OF  A 
CALL  LUDCMQ  (A,  MFP1,  10) 

C  RESET  THE  R.H.S.  INTO  C.  WE  NEED  TO  DO  THIS  FOR  EACH  DEGREE. 

MS PI  =  MS  +  1 
DO  95  I  =  MSP1 ,MFP1 
DO  90  J  =  1,1 
90  C(J)  =  A(J,MFP2) 

CALL  SOLNQ  (A,  C,  I,  10) 

IM1  =  1-1 

C  NOW  WRITE  OUT  THE  COEFFICIENTS  OF  THE  LEAST  SQ  POLYNOMIAL. 

WRITE(6 , 202)  IM1 ,  (C(J),  J  =  1,1) 

202  FORMAT (//,2X, 'FOR  DEGREE  OF  M2,1  COEFFICIENTS  ARE:  \/,10X, 

1  11F11.3) 

C  COMPUTE  AND  PRINT  THE  VALUE  OF  BETA  =  SUM  OF  DEVIATIONS  SQUARED, 

C  WHICH  IS  =  (N  -  1  -  M). 

BETA=0 . 0 

DO  94  IPT  =  1,N 

SUM=0 . 0 

DO  93  ICOEF  =2,1 

JCOEF=  I-ICOEF  +  2 
SUM=(SUM  +  C(JCOEF) )*X(IPT) 


CONTINUE 
SUM=SUM+C( 1 ) 

BETA=BETA  +  (Y(IPT)  -  SUM)**2 

CONTINUE 

BETA=BETA/(N-1) 

WRITE (6, 203)  BETA 

FORMAT (1H1,10X,  1  BETA  IS  ',F10.5) 

CONTINUE 

STOP 

END 


SUBROUTINE  LUDCMQ  (A,  N,  NDIM) 

DIMENSION  A(NDIM , NDIM) 

C  THIS  SUBROUTINE  FORMS  THE  LU  EQUIVALENT  OF  THE  SQUARE 
C  COEFFICIENT  MATRIX  A.  THE  LU,  IN  COMPACT  FORM,  IS  RETURNED 
C  IN  THE  A  MATRIX  SPACE.  THE  UPPER  TRIANGULAR  MATRIX  U  HAS  ONES 
C  ON  ITS  DIAGONAL  -  THESE  VALUES  ARE  NOT  INCLUDED  IN  THE  RESULT. 
DO  30  1=  1 ,N 
DO  30  J  =2 , N 
SUM=0 . 

IF(J.GT.I)  GO  TO  15 
JM1  =  J  -  1 
DO  10  K=1 , JM1 

10  SUM  =  SUM  +  A(I,K)*A(K, J) 

A(I , J)  =  A( I , J )  -  SUM 
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GO  TO  30 

15  IM1  =1-1 

IF(IMl.EQ.O)  GO  TO  25 
DO  20  K  =  1,1111 

20  SUtl  =  SUM  +  A(I,K)*A(I<,J) 

C  TEST  FOR  SHALL  VALUE  ON  DIAGONAL. 

25  IF(ABS(A(I, I)). LT. 0.0000001)  GO  TO  99 

A( I , J)  =  (A(I , J)  -  SUM) /A ( I , I ) 

30  CONTINUE 

RETURN 

99  WRITE ( 6 , 100 )  I 

100  FORMAT (/ , '  REDUCTION  NOT  COMPLETED  BECAUSE 
1  SMALL  VALUE  FOUND  FOR  DIVISOR  IN  ROW  ',13) 

RETURN 

END 


SUBROUTINE  SOLNQ  (A, B.N.NDIM) 
DIMENSION  A(NDI11,NDIM)  , B (NDIM) 


C  THIS  SUBROUTINE  FINDS  THE  SOLUTION  TO  A  SET  OF  N  LINEAR 
C  EQUATIONS  THAT  CORRESPONDS  TO  THE  RIGHT  HAND  SIDE  VECTOR  B 
C  THE  A  MATRIX  IS  THE  LU  DECOMPOSITION  EQUIVALENT  TO 
C  THE  COEFFICIENT  MATRIX  OF  THE  ORIGINAL  EQUATIONS,  AS 
C  PRODUCED  BY  LUDCMP.  THE  SOLUTION  VECTOR  IS 
C  RETURNED  IN  THE  B  VECTOR. 

C  DO  THE  REDUCTION  STEP 
B(1 )  =  B(l)/A(l,l) 

DO  20  I  =  2 ,N 
IM1  =  I  -  1 
SUM  =0.0 
DO  10  K  =  1 , IM1 
10  SUM  =  SUM  +  A(I,K)*B(K) 

20  B( I )  =  ( B ( I )  -  SUM)/A(I , I) 

C  NOW  WE  ARE  READY  FOR  BACK  SUBSTITUTION.  REMEBER  THAT  THE 
C  ELEMENTS  OF  U  ON  THE  DIAGONAL  ARE  ALL  ONES. 

DO  40  J  =  2 ,N 
MMJP2  =  N  -  J  +2 
NMJP1  =  N  -  J  +  1 
SUM  =  0 
DO  30  K 


DO  30  K  =  NMJP2 ,N 

SUM  =  SUM  +  A(NMJP1 ,K)*B(K) 

B(NMJPl)  =  B(NHJPl)  -  SUM 

RETURN 
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PROGRAM :  POLY 

AUTHOR:  BRUCE  F.  HILLARD 

DATE :  SEPTEMBER  1984 

PROGRAM  DESCRIPTION:  THIS  FORTRAN  PROGRAM  USES  DISSPLA 

SUBROUTINES  TO  GENERATE  PLOTS  OF 
MINI -RANGER  RANGE  ERRORS  THAT  ARE 
PREDICTED  FROM  POLYNOMIAL  REGRESSION 
CURVES.  ALSO  PLOTTED  ARE  THE 
MEAN  RANGES  ERRORS  USED  TO  CREATE 
THE  REGRESSION  CURVES.  THE  RANGE 
ERRORS  ARE  ASSUMED  TO  BE  A  FUNCTION 
OF  SIGNAL  STRENGTH. 

INPUT:  PLOT  VALUES  ARE  GENERATED  WITHIN  THE  PROGRAM  BY 

THE  REGRESSION  CURVES.  THE  MEAN  RANGE  ERRORS  ARE 
INPUT  THROUGH  A  DATA  STATEMENT. 

OUTPUT:  THE  RESULTANT  PLOT  CAN  BE  DIRECTED  TO  ONE  OF  THREE 

PLOTTERS:  VERSATEC  8222A,  VERSATEC  0900A,  OR  A 

TEKTRONIX  4631  HARDWARE  DEVICE. 


c************************************* ******************************** 


C 

IMPLICIT  REALM  (A-H.O-Z) 

DIMENSION  X(14) ,Y(14) ,Y1 (100) ,Y2(100) , Y3 (100) ,X1 (100 ) , 

1  X2(100) ,X3(100) , IPAK(IOO) 

DATA  X(1),X(2) ,  X(3) ,X(4) ,  X(  5) ,  X(6 ) ,  X(7) ,X(8) ,X(9) ,X(10) ,X(11) , 

1  X(12) ,X(13) , 

1  X(14)/4.0, 5. 0,6.0, 7 .0,8.0,9.0,10.0,11.0,12.0,13.0,14.0, 

1  15.0,16.0,17.0/ 

DATA  Y(l) ,  Y(  2 ) ,Y(3) ,Y(4) ,Y(5) ,Y(6) ,Y(7) ,Y(8) ,Y(9) , Y( 10 ) ,Y(11) , 

1  Y(12  ,Y(13), 

1  Y( 14) /2. 9, 0.0, -1.4, -1.3, -2. 6, -1.6, 0.0, 0.1, 0.1, 0.1, -0.1, 0.7, 

1  1.2, 1.1/ 

C 

C 

C  *************  PERFORM  THE  NECESSARY  PLOT/SETUP  ROUTINES. 

C 

C  **********  INITIALIZE  PLOTTER  AND  SET  PAGE  SIZE. 

C  CALL  PRTPLT(72 , 6 ) 

CALL  COMPRS 

C  CALL  VRSTEC  (0,0,0) 

C  CALL  TEK618 

CALL  PAGE  (11 . ,8.5) 

C  CALL  BLOWUP (1.25) 

C 

C  **********  MAKE  SS  AXIS  IN  INTEGERS 
C 


CALL  XINTAX 


C 

C  **********  ORIENT  Y  AXIS  NUMBERS  HORIZONTALLY. 

C 

CALL  YAXANG  (0.0) 

C 

c  **********  SELECT  TOTAL  PLOT  AREA. 

C 

CALL  AREA2D(3. ,6. ) 

C 

C  **********  SCALE  AND  LABEL  AXES  WITH  SELF-COUNTING  OPTION. 
C 

CALL  YNAME(  'RANGE  ERROR  (M)  $\100) 

CALL  XMAME  (’SIGNAL  STRENGTHS ‘ , 100 ) 

CALL  BASALF( ’STAND’ ) 

C 

C 

c  **********  DEFINE  BOTH  AXES  AND  FRAME  THE  SUBPLOT  AREA. 

C 


CALL  DUPLX 


CALL  GRAF (18.0,1.0,3.0,-4.0,1.0,10.0) 

CALL  THKFRM( . 02 ) 

CALL  FRAME 
C 

c************  CREATE  THE  LEGEND  THAT  IDENTIFIES  EACH  PLOT. 

C 

C 

C 

c  **********  OTHER  PLOT/LEGEND  INFORMATION: 

C 

C 

C 

c****  COMPUTE  THE  FIRST  DEGREE  VALUES. 

C 

9  Al=-0 . 883 

B1=0 . 079 
C 

TEMP 1=4.0 
20  CONTINUE 

DO  200  1=1,65 
XI (I )=TEMP1 
Y1(I)=  A1+(B1*X1 (I) ) 

TEMP1=X1 (I)+0 .20 
200  CONTINUE 
C 
C 

C****  COMPUTE  THE  SECOND  DEGREE  VALUES. 

C 

A2=4 . 540 
B2=-l . 133 
C2=0 .058 
C 

TEMP 2 =4 . 0 
30  CONTINUE 

DO  300  J=1 , 65 

X2(J)=TEMP2 

Y2  ( J ) =A2+ ( B2*X2 ( J ) ) + ( C2*X2 ( J ) **2 ) 

TEMP2=X2 ( J ) +0 . 20  ' 

300  CONTINUE 
C 
C 

C****  COMPUTE  THE  THIRD  DEGREE  VALUES. 

C 

A3=17 . 337 
B3=-5 . 661 
C3=0 . 530 
D3=-0.015 
C 

TEMP3=4 . 0 
40  CONTINUE 

DO  400  K=1 ,65 

X3(K  =TEMP3 

Y3(K)=A3+(B3*X3(K) )+(C3*X3(K)**2)+(D3*X3(K)**3) 
TEMP3=X3 (K)+0 . 20 
400  CONTINUE 
C 

c  ***********  SET  up  LEGEND  INFORMATION. 

c 

MAXLIN=LINEST ( IPAK , 100 , 25 ) 

1  CALL  LINESP (2.0) 

2  CALL  LINES ( 'COMPUTED  MEAN  VALUESS 1 , IPAK , 1 ) 

3  CALL  LINES ('FIRST  DEGREE  VALUESS ', IFAK , 2 ) 

4  CALL  LINESf 'SECOND  DEGREE  VALUESS ', IPAK , 3 ) 

5  CALL  LINES ( ' THIRD  DEGREE  VALUESS ', IPAK, 4) 

6  CALL  MYLEGN( 'CURVE  TYPES: ',12) 

7  CALL  LEGLIN 

C 

C 

C****  plot  THE  MEAN  ERRORS  AND  THE  COMPUTED  VALUES. 

C 


on  n  on  on  n  oon  on  n  noonn  no  n  n  o  ooooooooo  o  on  on  on 


CALL  THKCRV ( . 025 ) 
CALL  CURVE (X , Y , 14 , 0 ) 


CALL  THKCRV (.02 5) 

CALL  DOT 

CALL  CURVE (XI , Y1 , 65 , 0 ) 


CALL  THKCRV (.02 5) 

CALL  DASH 

CALL  CURVE (X2,Y2, 65,0) 


CALL  THKCRV (.02 5) 

CALL  CHNDOT 

CALL  CURVE (X3,Y3, 65,0) 

CALL  LEGEND ( I PAK , 4 , 1 . 0 , 4 . 0 ) 

**********  ALL  CURVES  HAVE  BEEN  PLOTTED. 
**********  NOW  COMPUTE  THE  M.A.S.S.  AND  PLOT  IT. 

CLAST=X3 ( 1 ) 

DO  500  K=2 , 65 

CHECK=ABS (X3(K) -CLAST ) 

IF ( CHECK. GE. 2.0)  GO  TO  600 
CLAST=X3 (K) 

500  CONTINUE 

********  NOW  PLOT  THE  M.A.S.S. 

600  CONTINUE 
L=K 

XF=X3 (L) 

XT=XF 

YF=Y3 (L)+0 . 5 
YT=Y3 (L)+0 . 1  • 

CALL  VECTOR ( XF ,YF , XT , YT , 1303) 

XX=XF-0 . 3 
YY=YF+0 . 7 

CALL  RLMESS ('M.A.S.S.S' ,100,XX,YY) 

CALL  MESSAG( 1  CODE  NO.  0$ ',100 , 5 . 0 , 5 . 0) 

CALL  MESSAG( ' CODE  NO.  1$ 1 , 100 , 5 . 0 , 5 . 0) 


**********  CLOSE  OUT  METAFILE  AND  RETURN  TO  LEVEL  ZERO. 

CALL  ENDPL(O) 

CALL  DONEPL 


STOP 
DEBUG 
AT  1 
TRACE 
END 


TRACE 

ON 
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PROGRAM:  BASELINE 
AUTHOR:  BRUCE  F.  HILLARD 

DATE :  AUGUST  1985 

PROGRAM  DESCRIPTION:  THIS  FORTRAN  PROGRAM  USES  DISSPLA 

SUBROUTINES  TO  GENERATE  PLOTS  OF 
MINI-RANGER  RANGE  ERROR  VERSUS 
SIGNAL  STRENGTH  BY  BASE  LINE. 

INPUT:  INPUT  COMES  FROM  A  DATA  FILE  REFERENCED  WITH  A 

FILEDEF. 

OUTPUT:  THE  RESULTANT  PLOT  CAN  BE  DIRECTED  TO  ONE  OF  THREE 

PLOTTERS:  VERSATEC  8222A,  VERSATEC  0900A,  OR  A 
TEKTRONIX  4631  HARDWARE  DEVICE. 

k  -k  -k  -k  •k'k  -k  -k  *  *  A  *  *  *  *  *  *  *  *  *  •k  *  *  *  *  *  *  *  *  *  "k  -k  *  *  A  *  -k  'k  *  ‘k  -k  -k  -k  -k  -k  x  -k  *  -k  A  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

INTEGER  M , A , B , C ,D , E , X, XX , I , P , Q , R , S , T ( SS , TIME , RAVE , CODE , ANT , BASEL 
1  SSA , SSB , SSC , SSD , SSE , ANTA , IPAK , ANTB , ANTC , ANTD , ANTE , SSX 
DIMENSION  BASEL( 1500 ) , 

1  ANT ( 1500 ) , CODE ( 1500 ) ,TIME ( 1500 ) , XRERRE (1500) , 

1  SS(1500) ,XRERR(1500) ,XRATE(1500) ,VRERR(1500) ,VRATE(1500) , 

1  XRERRX(1500) ,XRERRA( 1500) ,XRERRB( 1500) ,XRERRC( 1500) ,XRERRD( 1500) , 
1  SSA(1500) , SSB (1500) , SSC( 1500) , SSD ( 1 500 ) , SSE ( 1 500 ) , RAVE ( 1500 ) , 

1  TEMPA (1500), TEMPB ( 1 500 ) , TEMP C ( 1500), TEMPD (1500), TEMFE (1500), 

1  SSX( 1500 ), TEMPX(1500), T0TAL(1500) ,P( 1500), 0(1500 )  ,R( 1500) , 

1  T0TALA(1500)  ,TOTALB(  1500)  ,S(1500)  ,T(1500) , 

1  TOTALC ( 1500 ) ,TOTALD ( 1500 ) , TOTALE ( 1500 ) 

REAL  XRERR , XRERRA , XRERRB , XRERRC , XRERRD , XRERRE , XRATE , VRERR , 

1  TOTALA , TOTALB , TOTALC , TOTALD , TOTALE , TEMPX , 

1  TEMPA , TEMPB , TEMPC , TEMPD , TEMPE , TOTAL , XRERRX 

**********  READ  IN  THE  STAT  DATA  FROM  A  FILEDEF. 


1=0 

7  1=1+1 

READ(4,208)  SS ( I ) , XRATE ( I ) , VRATE (I ) , XRERR( I ) , VRERR( I ) , BASEL ( I )  , 
1 CODE ( I )  , ANT ( I ) , RAVE ( I ) , TIME ( I ) 

IF (SS ( I ) . EQ . 99)  GO  TO  50 
GO  TO  7 
50  CONTINUE 
N=I 


**********  AT  THIS  POINT,  ALL  THE  DATA  HAS  BEEN  STORED  IN  ARRAYS. 
**********  CALCULATE  TOTAL  MEAN  FOR  ALL  DATA/ CODE. 


XX=0 

DO  31  M=  4,18 

TOTAL (M)=0 
X=0 

DO  30  1=1, N 

IF (CODE ( I ) .ME . 0 )  GO  TO  30 
IF(SS(I) . NE . M )  GO  TO  30 
TOTAL  ( M ) =TOTAL ( M ) +XRERR ( I ) 

X=X+1 

30  CONTINUE 

XRERRX ( M ) =TOTAL (M) / FLOAT ( X ) 

SSX(H)=M 

XX=XX+X 

31  CONTINUE 

**********  DETERMINE,  BY  BASE  LINE,  THE  MEAN  RANGE  ERRORS. 


nnn  onononnnno  no 


DO  165  M=4 , 1 8 

P(M)=0 


2  H  =0 
R  M)=0 
S  M  =0 
T(I-1)=0 
T0TALA(M)=0 
T0TALB(M)=0 
T0TALC(M  =0 
T0TALD(M )=0 
TOTALE (M)=0 
DO  160  1=1, N 


160  CONTINUE 


IF (CODE ( I ) .NE . 0 )  GO  TO  160 
IF ( SS ( I ) . NE . M )  GO  TO  160 
IF ( BASEL ( I ) . NE . 1 )  GO  TO  60 
TOTALA(M)=TOTALA(M)+XRERR( I ) 
P(M)=P(M)+1 
GO  TO  160 

IF(BASEL(I).NE.2)  GO  TO  70 
T0TALB(M  =TOTALB(M)+XRERR( I ) 
Q(M)=Q(M)+1 
GO  TO  160 

IF ( BASEL ( I ) .NE . 3 )  GO  TO  80 
TOTALC(M)=TOTALC(M)+XRERR( I ) 
R(M)=R(M)+1 
GO  TO  160 

IF ( BASEL ( I ) .NE . 4 )  GO  TO  90 
TOTALD(M)=TOTALD(M)+XRERR( I ) 
S(M)=S(M)+1 
GO  TO  160 

TOTALE ( M ) =TOTALE ( M ) +XRERR ( I ) 
T(M)=T(M)+1 


(M 

=SS 

I) 

M, 

=SS 

I 

M 

=SS 

I 

M 

=SS 

I 

(M 

=SS 

I 

c**************  NOW  CALCULATE  THE  MEAN  RANGE  ERROR  FOR  EACH  SIGNAL 
STRENGTH  FOR  EACH  BASE  LINE. 

XRERRA ( M ) =TOT ALA (M) / FLOAT ( P ( M ) ) 

XRERRB ( M ) =TOTALB ( N }/ FLOAT ( Q ( M ) ) 

XRERRC (M ) =TOTALC ( M } / FLOAT  R ( M ) ) 

XRERRD l M ) =TOT ALD  H ) / FLOAT  S  M ) ) 

XRERRE (M)=TOTALE(M) /FLOAT (T(M) ) 

A=A+P(M) 

B=B+Q(M 
C=C+R  M 
D=D+S(M) 

E=E+T (M) 


165  CONTINUE 


*************  PERFORM  THE  NECESSARY  PLOT/SETUP  ROUTINES. 

**********  INITIALIZE  PLOTTER  AND  SET  PAGE  SIZE. 

CALL  PRTPLT(72 ,6) 

CALL  COMPRS 

CALL  VRSTEC  (0,0,0) 

CALL  TEK613 

CALL  PAGE  ( 11.  ,8.5) 

CALL  BLOWUP (1.25) 

**********  MAKE  SS  AXIS  IN  INTEGERS 


non  nnn  non  nnnn  non  non  nnno  non  non  non 


CALL  XINTAX 

**********  ORIENT  y  AXIS  NUMBERS  HORIZONTALLY. 
CALL  YAXANG  (0.0) 

**********  SELECT  TOTAL  PLOT  AREA. 


CALL  AREA2D(8. ,6. ) 

**********  SCALE  AND  LABEL  AXES  WITH  SELF-COUNTING  OPTION 

CALL  YNAME( 'RANGE  ERROR  (M)  $’,100) 

CALL  XNAflE  (’SIGNAL  STRENGTHS  MOO) 

CALL  BASALF( ’STAND’ ) 

**********  CREATE  THE  HEADING  AND  SET  THE  PRINT  TYPE. 

CALL  SWISSM 

CALL  HEADIN  (’MINI-RANGER  III  BASE  LINE  CALIBRATIONS’ 
**********  DEFINE  BOTH  AXES  AND  FRAME  THE  SUBPLOT  AREA. 
CALL  DUPLX 

CALL  GRAF (18.0,1.0,3.0,-6.0,2.0,10.0) 

CALL  THKFRM(.02) 

CALL  FRAME 

**********  SET  UP  LEGEND  DATA  AND  TEXT  FOR  LEGEND  NAME. 
CALL  LINESP  (2.0) 

CALL  LINES  ( ’ CUMMULATIVES 1 , IPAK , 1 ) 

CALL  LINES  (’A  -  1061$ 1 , IPAK , 2) 

CALL  LINES  (’B  -  2417$ ’, IPAK , 3 
CALL  LINES  (’C  -  4083$ 1 , IPAK , 4 
CALL  LIMES  (’D  -  7017$ 1 , IPAK , 5 
CALL  LINES  (’E  -  9861$ 1 , IPAK , 6 ) 


**********  FIND  SIZE  FOR  LEGEND. 

XW=XLEGND ( IPAK , 6 ) 

YW=YLEGND(lPAK,6) 

CALL  MYLEGN ( 1  BASE  LINE-LENGTH  (M)$',100) 
CALL  LEGLIN 

**********  pLoT  BASE  LINE  A  VALUES. 

CALL  MARKER (15) 

CALL  DOT 

CALL  THKCRV ( . 025 ) 

DO  195  M=4 , 18 

TEMPA ( M ) =FLOAT ( SS A ( M ) ) 

195  CONTINUE 

CALL  CURVE ( TEMPA, XRERRA, A, 1) 

***********  pL0T  BASE  LINE  B  VALUES. 

CALL  MARKER( 15 ) 

CALL  DASH 

CALL  THKCRV (.02 5) 

DO  196  M=4 , 18 

TEMPB ( M ) =FLOAT ( SSB (M ) ) 

196  CONTINUE 

CALL  CURVE ( TEMPB, XRERRB,B,1) 
************  plot  BASE  LINE  C  VALUES. 

CALL  MARKER(15) 
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100,1.5,1) 


non  nnnnn  nnnn  non  nnnn  non  non 


CALL  CHNDOT 
CALL  THKCRV( .025) 

DO  197  M=4 , 18 

TEMPC ( M ) =FL0AT ( SSC ( M ) ) 

197  CONTINUE 

CALL  CURVE ( TEMPC, XRERRC,C,1) 

*************  pL0T  BASE  T.INE  D  VALUES. 

CALL  MARKER(15) 

CALL  CHNDSH 
CALL  THKCRV( .025) 

DO  198  M=4 , 18 

TEMPD ( M ) =FL0AT ( SSD ( M ) ) 

198  CONTINUE 

CALL  CURVE ( TEMPD , XRERRD , D , 1 ) 
*************  plot  BASE  LINE  E  VALUES. 
CALL  MARKER (15) 

CALL  MRSCOD  2,6,5. ,2. ,2. ,7. ,2. ) 

CALL  THKCRV ( . 025 ) 

DO  199  M=4 , 18 

TEMPE (M)=FLOAT(SSE (M) ) 

199  CONTINUE 

CALL  CURVE ( TEMPE, XRERRE,E,1) 


************  pL0T  CUMMULATIVE  VALUES. 

CALL  MARKER(15) 

CALL  RESET  ( ' MRSCOD 1 ) 

CALL  THKCRV (.025) 

DO  200  M=4 , 1 8 

TEHPX(M)=FLOAT(SSX(M) ) 

200  CONTINUE 

CALL  CURVE ( TEMPX , XRERRX , XX , 1 ) 

**********  OTHER  PLOT/LEGEND  INFORMATION: 

CALL  MESSAG( ' CODE  NO.  0$ 1 , 100 , 4 . 5 , 5 . 0 ) 
CALL  MESSAGE  CODE  NO.  1$  1  , 100 , 4 . 5 , 5 . 0 ) 

**********  DRAW  THe  LEGEND. 

CALL  LEGEND ( IPAK ,6,1.0,3.75) 


**********  READ  FORMAT  STATEMENT. 

208  FORMAT (6X, 12 , 7X, F6 . 1 , 8X, F3 . 1 , 10X, F4 . 1 , 11X, F3 . 1 , 3X, II , II , II , 12 , 14) 

**********  CLOSE  OUT  METAFILE  AND  RETURN  TO  LEVEL  ZERO. 

CALL  ENDPL(O) 

CALL  DONEPL 
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